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Why China is important in advancing the field of 


primatology 
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ABSTRACT 


Over the past few decades, field studies conducted 
by Chinese  primatologists have contributed 
significant new theoretical and empirical insights into 
the behavior, ecology, biology, genetics, and 
conservation of lorises, macaques, langurs, snub- 
nosed monkeys, and gibbons. With the recent 
establishment and inaugural meeting of the China 
Primatological Society in 2017, China has emerged 
as a leading nation in primate research. Several 
research teams have conducted long-term studies 
despite the difficult challenges of habituating and 
observing wild primates inhabiting mountainous 
temperate forests, and the fact that some 80% of 
China’s 25-27 primate species are considered 
vulnerable, endangered, or critically endangered and 
are distributed in small isolated subpopulations. In 
going forward, it is recommended that primatologists 
in China increase their focus on seasonal 
differences in the social, ecological, physiological, 
and nutritional challenges primates face in exploiting 
high altitude and cold temperate forests. In addition, 
provisioning as a habitation tool should be minimized 
or eliminated, as it is difficult to control for its effects 
on group dynamics, patterns of habitat utilization, 
and feeding ecology. Finally in the next decade, 
Chinese primatologists should consider expanding 
the taxonomic diversity of species studied by 
conducting research in other parts of Asia, Africa, 
and the Neotropics. 


Keywords: China; Conservation; Primate research; 
Ecology 


INTRODUCTION 


Since the first study of a primate in the wild by Clarence 
Raymond Carpenter in 1931, research in field primatology has 
been conducted principally by scientists from the United States, 
Europe, and Japan; regions largely devoid of nonhuman 
primates (Sussman, 2010). Over the past 20 years this has 
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begun to change, with primate habitat countries such as Brazil, 
Mexico, Argentina, India, and China developing a cadre of 
highly-trained researchers and a national primatological society 
(http://www. internationalprimatologicalsociety.org/affiliatedsociet 
ies.cfm). Among this group, China stands out for its rapid 
economic development, the government’s commitment to 
funding high-quality scientific research, within-country access to 
state-of-the-art genetic, hormonal, imaging, geographic 
information system (GIS), and nutritional laboratories, as well 
as strong expertise in computer modeling and statistical 
analysis (King, 2004; Zhou & Leydesdorff, 2006). Moreover, 
China is home to an extremely diverse primate community that 
includes 25-27 species from four subfamilies (Lorisinae, 
Cercopithecinae, Colobinae, Hylobatinae) and 7 genera 
(Nycticebus, | Macaca, Rhinopithecus, Semnopithecus, 
Trachypithecus, Hoolock and Nomascus) (Estrada et al., 2017). 
These species differ markedly in body mass, social organization, 
mating systems, feeding ecology, and life history strategies. 
For example, group size among Chinese primates ranges from 
less than five to over 400 individuals (Qi et al., 2014). Similarly 
there are species described as solitary or inhabiting 
neighborhoods, whereas others live as socially monogamous 
pairs, one adult male — two adult female groups, harems or 
large one male multi-female breeding groups (OMUs), multimale- 
multifemale groups, and modular or multilevel societies 
composed of several OMUs plus an associated all male unit 
(Fan et al., 2010; Qi et al., 2014). Given this diversity, and the 
fact that studies of wild Chinese nonhuman primates have been 
conducted principally by in-country teams of scientists, Chinese 
primatologists are positioned to both empirically and theoretically 
advance the discipline of primatology in new and exciting ways. 


THE CHALLENGES OF CONDUCTING PRIMATE FIELD 
STUDIES IN CHINA 


Long-term studies vs. multigroup comparisons 
Despite the considerable challenges of conducting long-term 
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field studies, several research teams in China have collected 
data on the same primate group or population for a period of 
10-30 years (e.g., Tibetan macaques (Macaca thibetana), 
golden snub-nosed monkeys (Rhinopithecus roxellana), rhesus 
macaques (Macaca mulatta) , black crested gibbons 
(Nomascus concolor), black-and-white snub-nosed monkeys 
(Rhinopithecus biet, white-headed langurs (Presbytis 
leucocephalus)). This commitment has resulted in an extensive 
longitudinal database from which to distinguish and compare 
year-to-year variation in individual fertility, reproductive success, 
social interactions, and group dynamics with long-term patterns 
that may be indicative of important evolutionary processes. 
Other teams have studied multiple groups of the same species 
living in different habitats and can address critical research 
questions concerning which specific aspects of a species’ 
behavior and biology are conservative and remain constant 
across different environmental conditions and which traits are 
more plastic and vary (and to what degree they vary) under 
changing conditions (Zhao et al., 2011). The analysis and 
publication of these data will likely change current perspectives 
on many aspects of primate behavior, ecology, social systems, 
and mating strategies. 


Physiological challenges exploiting high altitude forests 

A number of Chinese primates inhabit high altitude temperate 
forests (including gibbons, leaf-monkeys, macaques, and 
snub-nosed monkeys) characterized by relatively long, cold 
winters. Although primates living in both tropical and 
temperate regions face significant ecological challenges in 
locating resources and consuming a nutritionally balanced 
diet, temperate-living populations encounter additional 
problems associated with hypoxia and the increased 
energetic costs (or reduced energy expenditure) of remaining 
thermoneutral during cold daytime and nighttime 
temperatures. The set of behavioral, developmental, genetic, 
nutritional, and physiological factors that enable nonhuman 
primates to successfully exploit temperate habitats are poorly 
understood, and Chinese primatologists are poised to take the 
lead in developing new theoretical frameworks to examine 
concepts of adaptation and adaptability in primate evolution. 


Ecological data 

Over the past decade, Chinese primatologists have 
increasingly published their research in high-impact journals 
(Fan & Ma, 2018). The majority of these studies have focused 
on primate social organization, social relationships, 
reproductive behavior, population genetics, and group size 
and composition. Only a limited number of studies have 
included a strong ecological component in the research. In 
part this relates to the fact that many primate species in China 
have extremely large home ranges (e.g., golden snub-nosed 
monkeys have home ranges of >10 km?) and presently are 
restricted in their distributions to mountainous regions and 
steeply sloped terrain unsuited to agricultural development. 
This makes the systematic collection of ecological data on 
tree species abundance, distribution, and resource 
productivity more difficult (unpublished data). In addition, 
several primate groups that have been the target of long-term 
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studies are semi-provisioned to facilitate habituation. This 
complicates the ability of researchers to test theories 
examining interrelationships between day range, patterns of 
habitat utilization, diet, and nutritional strategies and should 
be minimized or eliminated. 

Overcoming these challenges to promote greater emphasis 
on the feeding ecology and foraging strategies of Chinese 
primates is needed to better understand the causative effects 
of food availability on group size, feeding competition, and 
reproductive success. In the absence of detailed ecological 
information, researchers are constrained to offer explanations 
of their behavioral results that are based on untested 
assumptions present in the primate literature. For example, 
although it is tempting to assume that golden snub-nosed 
monkeys living in high altitude forests are food-limited during 
the long winter season, reduced food availability is not 
consistent with data indicating that golden snub-nosed 
monkeys live in large, cohesive, multilevel societies of several 
hundred individuals during all months of the year. Moreover 
the assumption of reduced food availability seems to 
contradict the fact that female nutritional requirements 
increase during pregnancy and female golden snub-nosed 
monkeys are pregnant during the winter. Finally, all group 
members must increase food intake during the winter to offset 
the additional energetic demands required to remain 
thermoneutral. Thus, either food available to golden snub- 
nosed monkeys is not limiting during the winter, or individuals 
of this species have evolved a set of successful dietary, 
digestive, or behavioral strategies to overcome this challenge. 

A stronger ecological focus would provide a constructive 
framework for identifying functional and adaptive relationships 
among primate behavior, anatomy, and ecology, as well as 
opening new avenues of inquiry that include species 
differences in spatial memory, decision-making, and cognition. 
Equally important, a stronger ecological focus will enable 
researchers to identify which species of trees and food items 
are most critical for conserving vulnerable primate 
populations. This would facilitate the success of projects 
involved in habitat restoration, the construction of forested 
corridors designed to facilitate migration and gene flow among 
isolated subpopulations, and habitat expansion into 
environmentally suitable and newly set aside protected areas. 


PRIMATE CONSERVATION 


At present some 80% of primate species in China are listed as 
vulnerable, endangered, or critically endangered by the IUCN 
(Estrada et al., 2017). We are at a historic moment in which 
China has the opportunity and the human capital of well-trained 
scientists to reverse decades of environmental degradation and 
biodiversity decline. Conservation must be a priority and a major 
component of all field and captive studies of Chinese primates. 


PRIMATOLOGY AS A SCIENTIFIC DISCIPLINE 


The rapid expansion of Primatology in China has been fueled by 
an exponential increase in the number of individuals’ receiving 


advanced degrees. Fan & Ma (2018) report that between 1984— 
2016 some 480 Chinese researchers received either a M.A. 
degree or Ph.D. in primatology. In recognition of this and other 
scientific achievements, in 2017 the International 
Primatological Society and the Chinese Government officially 
sanctioned the creation of the China Primatological Society 
(CPS). The CPS had its inaugural meeting in August 2017 in 
Xian, China. Some 200 scientists, directors of research 
institutes, students, and local governmental officials attended, 
including a small number of international primatologists from 
Australia, Canada, Japan, and the United States. Given the 
increased prominence of Chinese primatologists, the CPS is 
positioned to sponsor an official international, peer-reviewed, 
high-impact, scientific journal for the study of primate 
behavior, ecology, conservation, and evolution. 

It is clear that Chinese scientists have and will continue to 
play an increasingly important role in advancing the science of 
primatology in the coming decades. To facilitate this, | 
encourage research teams to expand their commitment to 
cross-species comparisons rather than each team focusing it 
efforts on only one or two taxa. Moreover, few research teams 
are currently studying the two extant loris species in China 
(Nycticebus bengalensis and N. pygmaeus), as well as several 
species of macaques (Macaca munzala, M. leucogenys, M. 
leonia, M. arctoides) and langurs (Semnopithecus schistaceus, 
Trachypithecus crepusculus, T. pileatus, and T. phayrei) (Fan & 
Ma 2018). Given the number of primatologists in China, field 
studies of these taxa, many of which are threatened, must be 
expanded. Finally, in order for Chinese primatologists to 
maximize their impact on the discipline, researchers should 
consider studying primates inhabiting other regions. There are 
some 504 species of primates worldwide (Estrada et al., 2017). 
Conducting research on species of primates from other parts 
of Asia, Africa, and the Neotropics will provide Chinese 
primatologists with the broadest and most informed 
comparative perspectives needed to rethink, refine, and test 
new theories in primatology. 
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ABSTRACT 


China is one of the most dynamic countries of 
the world and it shelters some amazing levels of 
biodiversity, including some very special primate 
species. However, primarily as a result of forest 
loss, most of which occurred in historical times, 
approximately 70% of China’s primate species have 
less than 3 000 individuals. Here | evaluate one 
road for future conservation/development that could 
produce very positive gains for China’s primates; 
namely forest restoration. | argue that for a large scale 
restoration project to be possible two conditions must 
be met; the right societal conditions must exist and 
the right knowledge must be in hand. This evaluation 
suggests that the restoration of native forest to support 
many of China’s primates holds great potential to 
advance conservation goals and to promote primate 
population recovery. 


Keywords: Conservation; Forest restoration; 
Regeneration; Primate population dynamics 


INTRODUCTION 


The world is changing rapidly and China represents one of 
the most dynamic countries on earth and it shelters some 
amazing levels of biodiversity (i.e., >30 000 species of vascular 
plants (behind only Brazil and Colombia), ~2 340 species of 
terrestrial vertebrates (Liu et al., 2003)). Globally, biodiversity 
is being lost at an accelerating rate, with current extinction 
rates approximately 1 000 times higher than background rates 
(Pimm et al., 2014). Recent estimates suggest that 11 000-58 
000 species are lost each year and that surviving vertebrate 
species have declined in abundance by 25% since 1970 
(Dirzo et al., 2014). Humans are clearly responsible for this 
accelerating loss of biodiversity, including the endangerment of 
primates. Between 2000 and 2012, 2.3 million km? of forest 
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was lost globally and in the tropics forest loss increased each 
year (Hansen et al., 2013). To put this in perspective, this 
area is approximately the size of Mexico. Global estimates of 
the extent of wildlife over-exploitation are very poor. However, 
Bennett et al. (2000) estimated that six million mammals were 
hunted annually in Malaysian Borneo. With respect to climate 
change, temperatures are predicted to increase by 1.5 °C by 
the end of the 21st century (IPCC, 2014) and using moderate 
greenhouse gas emission estimates, it is projected that by 
2100 75% of all tropical forests present in 2000 will experience 
temperatures that are higher than the temperatures presently 
supporting closed canopy forests (Peres et al., 2016; Wright et 
al., 2009). 

China follows some of these patterns, but many aspects 
how China has changed are unique; they represent different 
challenges and most importantly different opportunities. For 
example, since 2000, China’s Gross Domestic Product has 
increased by approximately 270% and China is now the world’s 
largest economy (Ahrends et al., 2017). The economy has to 
support a population approaching 1.4 billion people (Deng et 
al., 2015; Wei & Ye, 2014). Both the economic and population 
growth have come at an environmental cost. It is estimated that 
China has lost between 1.9 and 2.7 million km? of its original 
forest in the last 2 000 years (based on models of habitat 
suitability; Ahrends et al., 2017); this is an area approximately 
the size of the Democratic Republic of Congo. This has cost 
China in terms of primate diversity and population size. At 
least three species have been extirpated from China (Pygathrix 
nemaeus, Hylobates lar yunnanensis, Nomascus leucogenys) 
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and the Hainan (Nomascus hainanus) and Cao-vit gibbons 
(Nomascus nasutus) will almost certainly not see the turn of the 
century without very effective conservation action (Fan, 2017; 
Turvey et al., 2017). Gibbons were once found from Xi’an in 
central China east to Shanghai and all the way south to the 
border (Fan, 2017; Turvey et al., 2017; Zhou & Zhang, 2013); 


now they are isolated in a few forest fragments to the south. 


Even orangutans were found in southern China only 12 500 
years ago (Husson et al., 2009; Steiper, 2006), but they are no 
longer found on mainland Asia. A recent analysis considering 
22 of the 27 primate species in China (Fan & Ma, 2018) 
suggests that 15 of the species have less than 3 000 individuals 
and that 81% of the populations of all Chinese primates are 
declining (Estrada et al., 2017). These sorts of statistics clearly 
indicate that China has lost a great deal of ground in the battle 
to conserve its biodiversity, but they also illustrate that China 
has great potential in terms of primate conservation. 

The objective of this opinion article for this special issue on 
Primates and Primatology in China is to evaluate one road 
for future conservation/development that could produce very 


positive gains for China’s primates; namely forest restoration. 


There is clearly the need for restoration, as many of China’s 
primates are only found in small isolated forests; thus 
expanding their habitat and connecting fragments is clearly 
vitally needed. 


MAKING RESTORATION A VIABLE CONSERVATION 
STRATEGY 


For a large scale restoration project to be possible two 
conditions must be met. First, there must be the right 
societal conditions to make restoration possible and second 


the knowledge must be in hand to carry out such projects. 


In terms of the societal conditions, it appears that the 
timing is right for restoration projects. In November 1988, 
government of China enacted the "Law of the People’s 
Republic of China on the Protection of Wildlife" to facilitate 
the protection and management of wildlife. This law is 


the first truly comprehensive law to protect wildlife in China. 


However, since 1998, the Chinese government has enacted 
several national biodiversity conservation regulations, such 
as the Natural Forest Protection Project and Ecological 
Forest Compensation, which have been effective in improving 
environmental conditions in many areas (Ren et al., 2015; Xu et 
al., 2009). Government financing for protected areas has also 
increased following the launch of the Wildlife Conservation and 
Nature Reserve Construction Project and the Special Fund for 
Capacity Building of National-Level Nature Reserves. China is 
investing substantially in reforestation and tree planting efforts 
and this has totalled more than US$ 100 billion in the past 
decade alone (Ahrends et al., 2017; Li et al., 2013; Viña 
et al., 2016; Zhang et al., 2000; Zhang, 2015). China now 
has the world’s largest plantation area (approx. 800 000 km? 


(approximately the size of Mozambique, Ahrends et al., 2017). 


At the same time, China is trying to reduce pressures on natural 
forests through strict bans on logging in primary forests and a 
massive expansion of its forest reserves to a current total more 


than 2 500 reserves covering 1.6 million km? (this area includes 
17.1% of the country, which is approximately the size of Iran or 
twice the state of Texas). Lastly, over the last two decades 
there has been a large movement of people from rural areas 
(i.e., next to the reserves where China’s primates are found) to 
the cities. In fact, the urbanization rate rose from 17.9 to 52.6% 
between1978 and 2012 and currently more than half of China’s 
population live in cities (Deng et al., 2015). 

In terms of the right societal conditions to make restoration 
possible, an area still requiring a great deal of effort is 
that of hunting. In many areas of southern China, where 
forest cover is still substantial, primate populations can be 
dramatically reduced because of hunting (Harrison et al., 
2016). Even though some primates have Class | Protected 
Animals in the Chinese Wildlife Conservation Law and hunting 
guns have been outlawed and confiscated, illegal hunting still 
frequently occurs as it has been a traditional practise that is 
promoted by poverty in local communities, the use of wildlife 
for medicinal uses, and poor knowledge and enforcement of 
the laws (Fan et al., 2014). For example, in south-west 
Guangxi Province, Francois’ langur (Trachypithecus francoisi) 
populations declined by 90% between the early 1980s and 
early 2000s, at which time the total population size was 
estimated to be only approximately 300 individuals found in 
14 isolated populations. The researchers conducting the later 
survey concluded that the primary threat to the langur was 
hunting, primarily for traditional medicine (Li et al., 2007). 

The second requirement is that the knowledge must be 
available to effectively carry out a large restoration projects. 
Globally, there are only a handful of studies about the response 
of primate communities to forest regeneration; however, these 
studies suggest that forests, and the primate community they 
support, can rebound very rapidly when left to recover or 
encouraged to recover. For example, Baya & Storch (2010) 
surveyed a site in Korup National Park, Cameroon that was 
abandoned 7-8 years previously and found populations of all 
eight species of diurnal primates that occur in the region; in 
addition, sighting frequency in this recovering area was not 
significantly different from other sectors of the park (Linder, 
2008). In Kibale National Park, Uganda, seven years after an 
area of grassland was replanted with trees as part of a carbon 
offset program (Omeja et al., 2012), all species of diurnal 
primates were present in high numbers, including endangered 
red colobus and chimpanzee. Such studies give hope for the 
future. 

Within China a great deal of research has accumulated over 
the last two or so decades that provides exactly the type of 
information needed for restoration/conservation efforts. To start 
such efforts accurate information on the state of the species 
to be targeted must be known. There is extensive survey 
information on the current size and threats to the country’s 
primates (e.g., Chen et al., 2015; Cui et al., 2016). The 
synthesis of this information will be vital in determining the 
locations to be prioritized in restoration efforts; however, other 
local aspects, such as the willingness of the local population to 
participate in restoration and to not hunt primates must also 
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be considered. Next, information must be available on the 
habitat requirements of the species targeted for help from a 
forest restoration project and again, there are many detailed 
projects focusing on primate habitat use (Fan et al., 2009, 


2012; Guo et al., 2008; Li & Rogers, 2006; Liu et al., 2013). 


Of particular importance if the restoration project is to involve 
active replanting of trees is detailed information on the diet of 
the primates. With this information, either the food trees of 
the animals can be planted, or species with similar nutritional 
traits can be used in the restoration effort. There have been a 
number of high quality studies done on the nutritional ecology 
of China’s primates (e.g., Liu et al., 2013; Ma et al., 2017; Hou 
et al., 2018). To make restorations efforts more effective, it 
is critical to understand the severity and geographical nature 
of threats to primates in China other than habitat loss (e.g., 
live capture for trade, bushmeat, and tourism; Li et al., 2003; 
Xia et al., 2016; Yang et al., 2007; Zhu et al., 2013). In 
addition, information on conservation genetics is needed to 
determine the nature of corridors that can allow population 
mixing (Liu et al., 2015). With respect to conservation genetics, 
information concerning past interpopulational gene flow and 
landscape barriers on both short (e.g., satellite and aerial 
photography) and long time frames (e.g., river barriers), while 
rare (but see Wang et al., 2017), will be particularly valuable 
in determining the dispersal capabilities of primates relative to 
different types of barriers and can be used to provide guidance 
as to the nature of corridors that can be constructed to facilitate 
population mixing (Guo et al., 2010; Wang et al., 2015, 2017). 


WHAT REMAINS TO BE DONE FOR RESTORATION IN 
CHINA 


While this existing information is extremely valuable, what 
is lacking is data specific to restoration efforts and studies 
documenting behavioral patterns and responses of primates 
to the regeneration of forests. For restoration projects to be 
most effective, information on the survival strategies used in 
regenerating forests of all the major primate groups, including 


prosimians, macaques, colobines, and gibbons will be vital. 


While this is a formidable task, it also represents an exciting 
challenge to academics and Universities; one where the value 
of research can be illustrated to the government and public. 
What now remains to be done is to pull this societal 
potential and information together to facilitate large-scale 
forest restoration efforts that is critically needed for primate 
conservation. Since so many primate species in China are only 
hanging on as small remnant populations that are often only 
composed of a few groups (e.g., the cao vit gibbon population 
is estimated to only involve 18 groups occupying forest patches 
of only 2 000 hm? with only 3-4 groups in China (Fan et 
al., 2011)) the only way to effectively promote conservation of 
these primates in the wild is through restoration. This will likely 
require some sort of well financed coordinating agency that is 
able to rally national and international scholars, conservation 
organizations, government agencies, and the public to first 
provided the needed scientific information and then to use 
this information to promote both natural regeneration and 
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reforestation efforts on a very large scale. 
CONCLUSION 


Primates are very charismatic and are generally liked by both 
the Chinese people and international groups, thus they can act 
as “Flagship” or “Guardian Angel” species to promote the value 
of restoring these lands to a native forest (Bicca-Marques & 
De Freitas, 2010; Simberloff, 1998). It is my opinion that this 
is an exciting time to integrate restoration into conservation 
strategies to make informed and effective conservation and 
management decisions for the primates of China. 
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ABSTRACT 

China supports the richest non-human primate 
diversity in the northern hemisphere, providing an 
excellent opportunity for Chinese primatologists to take 
a leading role in advancing the study of primatology. 
Primatology in China began to flourish after 1979. To 
date, Chinese primatologists have published more 
than 1 000 papers in journals indexed by the Chinese 
Science Citation Database and the Web of Science 
Core Collection, and universities and academic 
institutions have trained 107 PhD students and 370 
Masters students between 1984 and 2016. In total, 
the National Science Foundation of China has funded 
129 primate projects (RMB 71.7 million) supporting 59 
researchers from 28 organizations. However, previous 
research has also shown obvious species bias. 
Rhinopithecus roxellana, Rhinopithecus bieti, and 
Macaca mulatta have received much greater research 
attention than other species. Researchers have also 
tended to continue to study the same species (55.2%) 
they studied during their PhD training. To promote the 
development of primatology in China, we suggest (1) 
the need for a comprehensive primatology textbook 
written in Chinese, (2) continued training of more 
PhD students, and (3) encouragement to study less 
well-known primate species. 


Keywords: Gibbon; Snub-nosed monkey; Leaf 
monkey; Macaque; Slow lories 


INTRODUCTION 


China supports the richest diversity of non-human primates 
(hereafter primate) in the northern hemisphere. More than 20 
species from three families currently reside in China (Table 1), 
with the discovery of new species and populations still adding to 
the list. Remarkably, two new species (Hoolock tianxing: Fan 


Science Press 


et al., 2017; Macaca leucogenys: Li et al., 2015) have been 
described and new populations of three species (Nomascus 
nasutus: Chan et al., 2008; Trachypithecus pileatus: Hu et al., 
2017; Rhinopithecus strykeri: Long et al., 2012) have been 
discovered in China since 2006. In total, 32 primate species 
taxa have been reported in China by different researchers 
(Table 1). The Chinese population of eastern hoolock gibbon 
(Hoolock leuconedys), once thought to live in western Yunnan 
(Fan et al., 2011), is now recognized as a new species, known 
as the skywalker hoolock gibbon (H. tianxing: Fan et al., 
2017). Consequently, there are currently no populations of 
eastern hoolock gibbon in China. Thus, excluding species with 
unconfirmed distribution (Pygathrix nemaeus, Pan et al., 2007) 
or disputed taxonomy (e.g., Nycticebus sp., Pan et al., 2007) 
and the two species occurring in Indian controlled areas of 
southeastern Tibet (Trachypithecus geei and Hoolock hoolock, 
Ji & Jiang, 2004; Jiang et al., 2017), there are 27 species of 
primates in China available for study (Table 1). 

Chinese primates live in diverse habitats and many 
populations represent the northern most distribution of the 
species, genus, or family. For example, A. bieti lives in 
dark coniferous forest above 3 000 m in Yunnan Province 
and Tibet; Trachypithecus leucocephalus survives in karst 
forests without surface water in Guangxi; Nomascus concolor 
ranges in mountain forests that can be covered by snow in 
winter in central Yunnan. In these harsh seasonal habitats, 
primates have evolved different strategies from their lowland 
conspecifics or close relatives to cope with ecological and 
social challenges, such as living in large groups. Snub-nosed 
monkeys (Rhinopithecus spp.) live in multi-level societies of up 
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to 400 individuals (Qi et al., 2014) and leaf monkeys 
(Trachypithecus crepusculus) living in seasonal montane 
forests on Mt. Wuliang form multi-male multi-female groups 
of more than 100 individuals, whereas their low-land relatives 
normally live in one-male groups of less than 30 individuals 
(Fan et al., 2015a). Crested gibbons (Nomascus nasutus, 
N. concolor, and Nomascus hainanus) live in groups with two 


Table 1 Number of primate species of China in different literature 


breeding females in seasonal forests in China, whereas other 
gibbon populations typically live in adult pair groups with only 
one breeding female (Fan et al., 2015b; Guan et al., 
2018). Populations living in colder or more temperate 
habitats provide excellent opportunities to study evolutionary 
and ecological differences in population genetics, nutritional 
ecology, endocrinology, cognition, and microbiome. 





This Protection 


wi 
N 
wo 
A 


Species 


English Chinese Masters PhD 
papers papers thesis dissertation grants (RMB 


No. of NSFC Total funding 





i oeg Mo na m (n) 10 000) 
Nycticebus sp. J 
Nycticebus pygmaeus J J J J J Class | 5 17 2 3 2 42 
Nycticebus bengalensis J J J J J Class | 1 11 1 1 0 0 
Nomascus nasutus J J J Class | 10 8 4 0 1 17 
Nomascus leucogenys J J J J J Class | 0 1 4 0 1 23 
Nomascus hainanus J J J J J Class | 14 11 3 4 2 113 
Nomascus concolor J J J J J Class | 28 13 10 6 4 155 
Hylobates lar J J J J J Class | 2 10 1 0 0 0 
Hoolock tianxing J J Class | 5 35 12 0 3 156 
Hoolock leuconedys J J J J 
Hoolock hoolock MA 
Trachypithecus shortridgei J J J J J Class | 2 4 1 0 0 
Trachypithecus pileatus J J Class | 0 0 0 0 
Trachypithecus phayrei J J J J J Class | 3 37 4 4 1 75 
Trachypithecus leucocephalus — ,/ J J J J Class | 30 39 20 8 8 353 
Trachypithecus geei J 
Trachypithecus francoisi J J J J J Class | 12 111 42 4 8 258 
Trachypithecus crepusculus J JV J Class | 0 0 0 0 0 0 
Semnopithecus schistaceus J J J J J Class | 0 9 1 1 0 0 
Rhinopithecus strykeri J JV J Class | 4 3 2 0 1 62 
Rhinopithecus roxellana J J J J J Class | 135 215 112 37 33 2051 
Rhinopithecus brelichi J J J J J Class | 15 25 T 0 3 85 
Rhinopithecus bieti J J J J J Class | 65 109 43 18 18 789 
Pygathrix nemaeus J 
Macaca thibetana J J J J J Class II 77 86 35 7 10 353 
Macaca munzala J J J Not assessed 0 0 0 0 0 
Macaca mulatta J JV J J J Class II 37 353 73 24 16 1 239 
Macaca leucogenys J JV J Not assessed 3 1 0 0 0 0 
Macaca leonina J J J J J Class | 12 14 6 0 0 0 
Macaca cyclopis J J J JV J Class | 35 8 29 5 0 0 
Macaca assamensis J J J J J Class | 6 70 11 1 1 62 
Macaca arctoides J J J J J Class II 4 59 4 0 0 0 
Total 22 23 26 29 27 





1: Ji & Jiang, 2004; 2: Pan et al., 2007; 3: Jiang et al., 2015; 4: Jiang et al., 2017 and this study. This table also shows the species bias in 
research. Total funding for each species was incomplete because it was calculated from the database of the National Science Foundation of 


China. The State Forestry Administration of China, Ministry of Science and Technology of China, and other organizations also support primate 


studies, but we did not have access to this information. 


At present, most primates in China are threatened by habitat 
loss and degradation, illegal hunting, and small population size 
(Fan et al., 2017; Liu et al., 2015). Except for M. mulatta, 
Macaca thibetana, Macaca cyclopis, and R. roxellana, all other 
species have total populations of less than 10 000 individuals 
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(Li et al., in preparation). Two gibbon species have recently 
become extinct in China (Nomascus leucogenys: Fan et al., 
2014; Hylobates lar: Grueter et al., 2009) and the Hainan 
gibbon has a total population of less than 30 individuals, being 
the most endangered primate species in the world (Bryant et al., 


2016). According to the Chinese Wildlife Conservation Law, 22 
species are listed as National Class | Protected Species, three 
species (M. thibetana, M. mulatta, and Macaca arctoides) are 
Class II, and two species (Macaca munzala and M. leucogenys) 
have not yet been assessed (Table 1). Conservation oriented 
research can provide an informed scientific foundation to 
support species protection and preservation. 

Primatology is the scientific study of both living and extinct 
primates in either their natural habitats by field surveys or in 
laboratory experiments to understand aspects of their evolution 
and behavior (https://en.wikipedia.org/wiki/Primatology). Ji & 
Jiang (2004) dated the rise of primatology in China back to 
1862 when Swinhoe published his paper about mammals from 
the island of Formasosa (Swinhoe, 1862); however, Chinese 
scientists first studied primates in the 1950s (Shou, 1957; Tan, 
1957). Xia & Zhang (1995) and Ji & Jiang (2004) provided 
good summaries on the research and conservation of primates 
in China. Primatology, especially studies of wild populations 


in China, has experienced a strong development trajectory. 


Several teams have established field stations and conducted 
continuous long-term research for more than 20 years. 

The Chinese Primatological Society was formally 
established in 2017. To celebrate this landmark event, we 
organized this special issue on “Primates and Primatology 
in China”, which includes studies on gibbons, leaf 
monkeys, snub-nosed monkeys, and macaques, as well as 


commentaries from Prof. Paul A. Garber and Prof. Colin A. 


Chapman. To provide a basic introduction to primates and 
primatology in China, we briefly review the development of 
primatology in China through three indices: that is, number 
of publications, student training, and funding received from 


the National Natural Science Foundation of China (NSFC). 


We also consider the current challenges facing primatology 
development in China. This review excludes studies on extinct 
species or on primates as medical models. 


MATERIALS AND METHODS 


We searched for papers published in Chinese journals 
using the Chinese name of each extant species (including 
title, keyword, and abstract) within the Chinese Science 
Citation Database (CSCD) using Chaoxing, a Chinese 
literature database (http:/Awww.chaoxing.com). We searched 
for papers published in international journals through the Web 


of Science Core Collection (WSCC) database (http://apps. 


webofknowledge.com) using the Latin name of each extant 
species and refined the results by confining the county/region 
to “China”, “Hong Kong”, or “Taiwan”. We searched for Masters 
theses and PhD dissertations in Chinese within several 


databases, including CNKI, Wanfang (http:/Awww.wanfangdata. 


com.cn/), Airitilibrary (http://www.airitilibrary.cn), and Chaoxing 
(http:/Awww.chaoxing.com/), using the Chinese name of each 
species as the search term. We searched for grants for 
each species funded by the NSFC since its establishment 


in 1986 (http:/Awww.medsci.cn/sci/nsfc.do, http://nsfc.biomart. 


cn/). Although primate research has received funding from 
the State Forestry Administration of China, Ministry of Science 


and Technology of China, and other organizations, the NSFC 
is the largest funding source for research and the only open 
source to which we had access. We also used genus 
and family names when searching for papers published in 
Chinese and English, and for grants. We checked each 
paper and thesis and excluded those using primates as 
medical models. We also excluded grants for studies of 
more than one species (n=12) when determining total NSFC 
funding for each species. We calculated the proportion of 
international cooperation by dividing the number of papers 
in which an international organization was involved based 
on author affiliation by the total number of papers published. 
We also obtained anecdotal information from primatologists in 
the WeChat group “China Primatological Society” by asking 
“who teaches primatology-related courses to undergraduate 
students” and “who obtained a PhD degree aboard”. We 
analyzed data and constructed figures using Microsoft Excel 
2010. 


RESULTS 


Number of publications 

In total, we found 496 papers published by Chinese authors 
from the WSCC and 999 papers from the CSCD (http://www. 
chaoxing.com). Only 15 papers were published in CSCD 
journals and three papers in WSCC journals before 1979, 
after which primate studies began to flourish. We observed 
stable development of Chinese primatology between 1979 and 
2003, with China then entering a new phase after the 2002 
International Primatological Society meeting held in Beijing. 
After 2003, Chinese primatologists published a greater number 
of papers in international journals (WSCC database), and the 
number of papers in the WSCC database exceeded the number 
of papers published in the CSCD for the first time between 2014 
and 2017 (Figure 1). 

Chinese primatologists have also actively engaged in 
international cooperation since the rise of primatology in 
China. In total, scholars based at international universities and 
academic institutions have been involved in 228 English papers 
(46.0%, n=496). The proportion of international cooperation 
in WSCC papers has remained above 25% over time (Figure 
1). In contrast, very few (3.3%, n=999) CSCD papers have 
included international scholars (Figure 1). 


Student training 

China commenced student training in primate research after 
1979. However, primatology is not a distinct discipline in China, 
with graduate students who have studied primates receiving 
their degrees in zoology, ecology, genetics, physiology, and 
anthropology. The first graduate student (Ya-Wen Huang) 
who studied the habitat of the Hainan gibbon obtained her 
Masters degree at Sun Yat-Sen University in 1984. The first 
PhD student (Ya-Ping Zhang) studied genetic diversity of the 
genus Macaca and obtained his PhD degree from the Kunming 
Institute of Zoology, Chinese Academy of Sciences (CAS) in 
1991. Since then, there has been a dramatic increase in 
the number of graduate students studying primates (Figure 
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2). From 1979 to 2016, 109 graduate students have obtained 
a PhD and 370 students have obtained their Masters in 
primate-related research. From 2013 to 2016, an average of 
eight PhD students and 31 Masters students have received 
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their degree each year (Figure 2). Based on information 


collected from the WeChat group “China Primatological 
Society”, six Chinese primatologists have obtained their PhD 
degree abroad. 
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Figure 1 Number of papers published per year by Chinese authors 


From the Web of Science Core Collection (WSCC: http://www.webofknowledge.com) and Chinese Science Citation Database (CSCD: http://www.chaoxing.com), 


and proportion of international cooperation during each time period (The 2017 database is not complete). 


The number of organizations that train Masters and PhD 
students in primate studies has also increased in China (Figure 
2). Prior to 2000, only four organizations (Kunming Institute 
of Zoology (CAS), Peking University, Beijing Normal University, 
and Institute of Zoology (CAS)) trained PhD students and seven 
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organizations trained Masters students in the study of primates. 
These numbers increased to 33 and 63, respectively, from 2014 
to 2016 (Figure 2). Currently, however, only two professors 
teach primatology to undergraduate students in China. 





1979-83 1984-88 1989-93 1994-98 1999-03 2004-08 2009-13 2014-16 
Time period (year) 


== Master training organization =E= Doctoral training organization 


=< Master thesis per year 


= PhD dissertation per year 


Figure 2 Average number of students per year that studied primates as part of their research, and the accumulated number of 


training organizations in China 
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NSFC funding 

Although the NSFC was established in 1986, it only funded 
its first primate field project in 1996 (study commenced in 
January 1997). Further studies were funded after 2009 
(Figure 3). In total, the NSFC has funded 129 primate 
projects (RMB 71.7 million), supporting 59 researchers from 
28 organizations. Of the 129 projects, 112 focused on a 
single species, 12 on a single genus, two at the family level, 
and three on the order Primates. Ecology and conservation, 
social behavior, and genetics and genomics are the three 
most common research fields, with several projects studying 
cognition, microbiome, anatomy, and physiology, and three 
focusing on education, computer facial recognition, and captive 
management (Figure 4). 
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Figure 3 Number of grants related to primate research funded 
by the National Natural Science Foundation of China 
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Figure 4 Research fields of the 129 NSFC-funded primate 
projects 


Species bias 

We found an obvious species bias in primate research in 
China (Table 1). Two endemic snub-nosed monkey species 
(R. roxellana and R. bieti) and common Macaca mulatta have 
received more research attention in China than any other 
taxa. More students have studied these species as their 
thesis subject, more projects on these species have been 
funded, and more papers have been published. At the 
same time, other species such as Trachypithecus crepusculus, 
Trachypithecus pileatus, Semnopithecus schistaceus, and 


Nycticebus bengalensis have been completely overlooked 
(Table 1). 

Thirty-three researchers who studied primates for their PhD 
received NSFC grants. However, except for a few senior 
professors, only two Masters degree holders have received 
grants from the NSFC. Researchers with a PhD have tended to 
continue studying the same species (55.9%), genus (70.6%), 
or family (85.3%). 


DISCUSSION 


Chinese scientists have made significant progress in primate 
research and student training since 1979. However, both total 
funding and total number of PhD students remain small. China 
has only trained 109 PhD students since 1979, with less than 
eight PhD students defending their degree each year from 
2014 to 2016. As Masters students have rarely been funded 
by the NSFC, and therefore have been unable to develop 
a successful research career, we must train and encourage 
more PhD students in order to promote primate research and 
conservation in China. 

We also need to pay attention to the obvious species bias 
shown in current research. Some species like R. roxellana 
and R. bieti are well studied and their populations are stable or 
increasing, whereas species like T. pileatus and N. bengalensis 
have received little attention and are on the edge of extinction. 
A pattern of continuing to study the same species after 
PhD training has only exacerbated the situation. Therefore, 
researchers should be encouraged to study overlooked species 
or focus on inter-species comparisons, as most NSFC projects 
have focused on single species. 

Another issue facing Chinese students is the lack of Chinese 
language textbooks on primatology. Currently, students must 
read English texts when learning and developing an interest 
in primatology, which can inhibit and discourage those with 
poorer English skills. Although many professors teach zoology 
or ecology in China, at present only Dr. Peng Zhang 
teaches primatology at Sun Yat-Sen University, and Drs. 
Dong-Po Xia and Bing-Hua Sun teach primatology at Anhui 
University. Consequently, systematic training in primatology 
for undergraduate and graduate students remains deficient 
in China, which can increase the difficulties for post-doctoral 
researchers to expand their area or species of interest. Given 
that developing a new discipline in China is challenging, 
we urgently need a Chinese language textbook that will 
encourage interested students to learn about primatology as 
undergraduates. 

China has made significant progress in primatology 
since 1979. To promote further development, however, 
we need to establish a high-quality Chinese language 
primatology textbook, train more PhD students, and encourage 
post-doctoral researchers to study less well-known primate 
species in China. Furthermore, new students need to be 
inspired to study primates and public awareness of primatology 
and primate conservation should be increased. Chinese 
primatologists need to be encouraged to write popular science 
books and introduce primates during lectures and talks. We 
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sincerely hope that primates, primatology, and primatologists 
will have a brighter future in China. 
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ABSTRACT 


Gibbons in China represent the northernmost margin 
of present day gibbon species distribution (around 
N25°). Compared to tropical habitats, northern gibbon 
habitats are characterized by low temperatures and 


remarkable seasonal variation in fruit abundance. 


How gibbons adapt to their cold and seasonal 
habitats and what ecological factors affect their 
sociality are key questions for understanding their 
ecology and social system evolution, the elucidation 
of which will contribute to the conservation of these 
special populations/species. According to preliminary 
short-term studies, northern gibbons consume more 
leaves and use larger home ranges than tropical 
gibbons. Interestingly, some Nomascus groups 
consist of more than one adult female. However, 
these preliminary results are not well understood or 
incorporated into current socio-ecological theories 
regarding gibbon species. To better understand 
northern gibbons, our team has systematically studied 
three habituated groups of Nomascus concolor, three 
groups of N. nasutus, and two habituated groups 
of Hoolock tianxing since 2002. In this paper, 
we stress the challenges facing gibbons living in 
northern habitats and summarize their behavioral 
adaptations to their harsh environments. We also 
describe the northern gibbon social system and 
discuss the potential relationships between their 
ecology and sociality. Finally, we highlight future 
research questions related to northern gibbons in 


Science Press 


China. 


Keywords: Northern gibbon; Ecology; Social system; 
Adaptation; Group size 


INTRODUCTION 


Gibbons (Hylobatidae) are small apes that live in tropical 
and subtropical forests in southeastern Asia (Bartlett, 2011). 
These apes are characterized by several interrelated features, 
including long arms for brachiation, terminal branch feeding, 
frugivory, pair-living, and territoriality (Brockelman, 2009). In 
tropical forests, gibbons rely heavily on ripe fruit pulp (Chivers, 
1984), with fruits including figs comprising >60% of their annual 
diet, though siamang (Symphalangus syndactylus) is more 
folivorous (reviewed in Bartlett, 2011). Fruit is often distributed 
in small and defensible patches, resulting in food competition 
among females. Female gibbons are considered intolerant to 
each other and live solitarily (Bartlett, 2011; Kappeler et al., 
2013; Leighton, 1987; Wrangham, 1987), with males unable 
to successfully defend ranges occupied by two or more females 
(Bartlett, 2009). Males may thus benefit from pair-living through 
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1) infanticide defense against non-resident, unmated males 
(Borries et al., 2011; Opie et al., 2013; van Schaik & Kappeler, 
2003); 2) increased paternal certainty by mate guarding 
(Palombit, 1996, 1999); 3) limited opportunities for the female 
to evaluate potential sexual partners (Palombit, 1996, 1999); 
and 4) enhanced foraging efficiency by sharing knowledge with 
females (Brockelman, 2009) or by ensuring that resources are 
not unknowingly depleted by foraging apart (Terborgh, 1983). 
Females benefit from pair-living through food and territory 
defense, direct infant care provided by males (reported in 
siamangs, Lappan, 2008), and enhanced foraging efficiency 
(Brockelman, 2009). The mated pair use a small home 
range of on average 42 hm?, and actively defend 82% of 
the area (Bartlett, 2011) through direct intergroup conflicts 
(Bartlett, 2003) and regularly occurring loud morning songs 
(Geissmann, 2002; Mitani, 1985). Both sexes benefit from 
knowledge of this intensively-used territory and goal-directed 
foraging strategies (Brockelman, 2009). This, in turn, reinforces 
territory defense and small home range, and consequently 
often limits gibbon group size (Brockelman, 2009; Brockelman 
et al., 2014). Bartlett (2009) first integrated seasonality into 
territory defense. He argued that stable food abundance is a 
prerequisite for territory defense. However, most information 
on gibbons comes from only a few studied species (mainly 
Hylobates lar and S. syndactylus) living in lowland tropical 
forests (Bartlett, 2009, 2011; Brockelman et al., 1998; Chivers, 
1974; Lappan, 2008; Palombit, 1996, 1999). 

Historically, gibbons were distributed widely in China. As 
recently as 500 years ago, their melodic songs could be heard 
along The Three Gorges in central China and even further north 
(close to N35°) (Turvey et al., 2015; Zhou & Zhang, 2013). 
Although the known gibbon range has shrunk considerably 
from north to south, and is now restricted to Yunnan, Guangxi, 
and Hainan in southern China (Fan, 2017; Turvey et al., 
2015; Zhou & Zhang, 2013), these sites still represent the 
northernmost margin of present day gibbon distribution (around 
N25°, e.g., N. concolor in Wuliang Mountain: Jiang et al., 2006; 
H. tianxing in Gaoligongshan: Fan et al., 2011b, 2017), second 
only to the Hoolock gibbons in India and Myanmar (Groves, 
1967). In this paper, we arbitrarily refer to the populations of 
N. concolor, N. leucogenys, N. hainanus, N. nasutus, and H. 
tianxing in China as “northern” gibbons (Figure 1). It does not 
mean that these gibbons are distributed in northern areas but 
is used for comparison between these gibbon populations and 
gibbons living in “southern” tropical forests. 

Based on short-term studies of unhabituated northern 
gibbons in China, early researchers reported on the differences 
in ecology and behavior compared to tropical gibbons. 
For example, northern gibbons were less frugivorous (N. 
leucogenys: Hu et al., 1990; N. concolor: Chen, 1995; Lan, 
1993), occupied larger home ranges (N. hainanus: Liu et al., 
1989; N. leucogenys: Hu et al., 1990; N. concolor: Jiang et 
al., 1994a; Jiang & Wang, 1999), and some groups included 
more than one adult female (N. concolor: Haimoff et al., 1986, 
1987; N. hainanus: Liu et al., 1989; Jiang et al., 1994a, 
1994b; Jiang & Wang, 1999). Although early studies reported 
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obvious differences, detailed observations on these gibbons 
were hampered in China by a lack of habituated wild gibbons 
due to intense human hunting pressure and difficult terrain 
(Yin et al., 2016). Systematic, longitudinal research became 
possible only recently after several N. concolor and H. tianxing 
groups became habituated to observers in Mt. Wuliang and 
Mt. Gaoligong (Fan & Jiang, 2008; Fan et al., 2008, 2009a, 
2013; Yin et al., 2016). In 2006, N. nasutus was rediscovered 
in a small karst forest patch in Bangliang, Jingxi County, 
Guangxi (Chan et al., 2008). Researchers could observe these 
gibbons from fixed observation posts located on hill tops or 
forest gaps due to the steep karst terrain and low canopy 
height (Fan et al., 2011a, 2012). With accumulating data 
from different species living in different habitats, the picture of 
northern gibbon ecology and sociality, although still incomplete, 
is becoming clearer. 


CHALLENGES FACING NORTHERN GIBBONS IN CHINA 


The human population in China has increased rapidly and 
currently stands at over 1.3 billion (National Bureau of 
Statistics of China, http://www.stats.gov.cn/english/). The 
resulting increase in demand for cropland, together with the 
steel-making campaign of 1958, has led to the clearance 
of vast tracts of forest (Jiang et al., 2006; Zhang, 2000). 
Deforestation and habitat fragmentation caused by expansion 
of crop cultivation and plantations are the main causes of 
gibbon decline and extinction in China (Fan et al., 2011a, 
2011b, 2013, 2014; Fan, 2017; Jiang et al., 2006; Zhang et 
al., 2010; Zhou et al., 2005), with large-scale hunting in the 
1970s and 1980s also having a considerable impact. A total 
of 81 365 primate skins were collected by Chinese government 
officials between 1974 and 1980 (Quan et al., 1981). While 
many skins were from macaque species (Macaca mulatta and 
M. arctoides), gibbons were also hunted during this period 
(Fan, 2017; Quan et al., 1981). In the late 1980s and early 
1990s, many Chinese gibbon populations had declined to such 
low levels that recovery was considered unlikely. By the 
late 1980s, the remaining populations of three species had 
decreased to less than 40 individuals (N. hainanus: Liu et 
al., 1989; H. lar: Yang et al., 1987; N. leucogenys: Hu et al., 
1990), and until recently N. nasutus was considered extinct 
in China (Tan, 1985). To date, the numbers of H. lar and 
N. leucogenys have not recovered and both species are now 
assumed to be ecologically extinct in China (Fan et al., 2014; 
Grueter et al., 2009). After more than 20 years of intensive 
conservation efforts, N. hainanus remains highly threatened 
with extinction (Fellowes et al., 2008; Zhou et al., 2005) and is 
the most endangered primate species in the world (Fan, 2017), 
with only one population of less than 30 individuals surviving 
in Bawangling. Fortunately, a small population of about 20 
individuals of N. nasutus was rediscovered in Bangliang, Jingxi 
County, Guangxi in 2006 (Chan et al., 2008). The remaining 
gibbon populations in China survive in montane or karst forests 
in several isolated mountain regions that are unsuitable for 
cultivation. 

Gibbon habitats in these northern mountain areas are much 


colder than tropical forests (Table 1). For example, the annual 
mean temperature in Dazhaizi, Mt. Wuliang is 16.6 °C (at 
2 000 m as.l.: Fan & Jiang, 2008) and is even colder at 
Nankang (13.3 °C, at 2 200 m a.s.l.: Fan et al., 2013) and 
Banchang, Mt. Gaoligong (13.0 °C, at 2 300 m a.s.l.: Fan et al., 
unpublished data). At these sites, temperatures can drop below 
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zero and forests are sometimes covered in snow (Fan & Jiang, 
2008; Fan et al., 2013). Thus, the resident gibbons likely face 
increased thermoregulatory costs compared to tropical gibbon 
populations. Nakayama et al. (1971) reported that the energy 
expenditure of outdoor-living captive Japanese macaques (M. 
fuscata) at 5.2 °C was 2.5 times greater than that at 29.5 °C. 


Figure 1 Distribution of the gibbon species: Nomascus (7 species); Hoolock (3 species); Hylobates (9 species); Symphalangus (1 
species) (According to Thinh et al., 2010a, 2010b; Fan et al., 2017; IUCN, 2017) 


Our gibbon research sites are also marked on the map. 


Table 1 Behavioral ecology information of northern gibbons 
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Fruit productivity decreases with increasing latitude and 
altitude (Hanya & Aiba, 2010), and northern gibbons are known 
to consume less fruit than tropical gibbons (reviewed in Fan et 
al., 2013). Fruiting seasonality also increases with latitude, and 
fruiting seasons are longer in the tropics and shorter toward the 
poles (Ting et al., 2008). At Mt. Ailao, a site supporting one of 
the two largest populations of N. concolor in China (Li et al., 
2011), figs are absent and other fruit is only available for four 
months of the year (Chen, 1995; Sun et al., 2007). Haimoff 
et al. (1987) argued that Mt. Wuliang represents an extreme 
environment for gibbons due to the unavailability of fruit for 
much of the year, during which time gibbons must subsist on 
leaves (also see Fan et al., 2009a). In Mt. Gaoligong, Nankang, 
fruit availability drops close to zero between December and 
February (Fan et al., 2013). These findings raise the question 
of how northern gibbons adapt to their cold and seasonal 
habitats. 


ECOLOGY OF NORTHERN GIBBONS 


Northern gibbons live in sizable home ranges (N. concolor: 
>100 hm?, N. nasutus: 130 hm?, H. tianxing: >93 hm? in 
Table 1; N. leucogenys: 540 hm?, Hu et al., 1990; N. hainanus: 
200-500 hm?, Liu et al., 1989, yearly 149 hm? in Bryant et 
al., 2016), which are considerably larger than those of tropical 
gibbons (reviewed in Bartlett, 2011). In mountainous areas, 
trees are patchily distributed in different microhabitats due 
to variables such as soil, altitude, sunlight, water flow, and 
topography (Fan et al., 2011a; Sun et al., 2007; Tian et al., 
2006). Tree species such as Saurauia napaulensis and Ficus 
neriifolia occur along streams on damp shady slopes, whereas 
Tetrastigma delavagi occurs on the sunny slopes at Dazhaizi, 
Mt. Wuliang (Tian et al., 2006). Ficus glaberrima is distributed 
on slopes and valleys, whereas Platycarya longipes occurs on 
hill tops. The karst habitat used by N. nasutus was classified 
into four different types (Fan et al., 20114). 

Patchily distributed plant species create variable food maps 
for northern gibbons. To obtain sufficient sustenance, the 
group of N. concolor in Mt. Wuliang (Fan & Jiang, 2008) 
and H. tianxing in Mt. Gaoligong (Zhang et al., 2014) 
shift their monthly home range according to the seasonal 
availability of food species. Although their monthly home 
ranges (N. concolor: 37+12 hm?; H. tianxing: 28+15 hm?) are 
comparable to the home ranges of tropical gibbons (avg. 34-49 
hm?, calculated from Table 17.3 in Bartlett, 2011, excluding 
Nomascus), they occupy a larger total home range to satisfy 
their annual food requirements (Fan & Jiang, 2008; Zhang et 
al., 2014). Furthermore, they do not use their large home 
range evenly, but instead heavily utilize smaller areas with 
greater food resources (Fan & Jiang, 2008). They also do 
not significantly increase their daily path length compared 
with tropical gibbons (northern gibbons in Table 1; tropical 
Hylobates gibbons: avg. 1 300 m), except for some siamang 
species (e.g., S. syndactylus: 0.8 km, reviewed in Bartlett, 
2011). 

In response to the overall low availability of fruit, early 
researchers suggested that northern gibbons ate more leaves 
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and less fruit than tropical gibbons (Chen, 1995; Hu et al., 
1990; Lan, 1993). This hypothesis was supported by several 
later studies with at least one-year duration. For example, 
the overall diet of N. concolor at Dazhaizi, Mt. Wuliang, was 
comprised of 46.5% leaves and 44.1% fruit (Fan et al., 2009a), 
with similar results reported for H. tianxing in Nankang, Mt. 
Gaoligong (49.1% fruit and 43.3% leaves, Fan et al., 2013) 
and N. nasutus in Bangliang (50% fruit and 35% leaves, Fan 
et al., 2012). Comparatively, the diets of tropical gibbons are 
typically comprised of 64% fruit and 25% leaves, except for the 
more folivorous siamangs (49% fruit and 38% leaves) (Bartlett, 
2011). Although northern gibbons consume less fruit than 
tropical gibbons, they prefer fruit when it is available (Fan et al., 
2012). When non-fig fruit is more abundant, H. tianxing feed 
on fewer leaves (Fan et al., 2013). Although Fan et al. (2009a) 
did not monitor fruit availability at Dazhaizi, Mt. Wuliang, the 
N. concolor group studied spent more time feeding on fruit in 
March of the second year when fruit was more available (Fan 
et al., 2009a). When fruit is less accessible or absent, northern 
gibbons switch their diet to leaves and buds (Fan et al., 2009a, 
2012, 2013), accounting for 83% of the N. concolor diet in 
March (Fan et al., 2009a), and increase their feeding time (Fan 
et al., 2008). 

When leaves dominant their diets, gibbons decrease their 
travel time and shorten their daily path length as leaves contain 
fewer available calories per gram than fruit (Raemaekers, 1978) 
and are more abundant and evenly distributed in the forest (Fan 
et al., 2008). This travel strategy has also been adopted by 
tropical gibbons (Chivers, 1974; Bartlett, 2009; Raemaekers, 
1978), with two groups of H. lar found to travel on average 1 
330 m (range from 620 m to 1 930 m, n=12 months) and 1 160 
m (range from 720 m to 1 660 m, n=12 months), they traveled 
short distance when dominantly feeding on leaves (Bartlett, 
2009). However, variation in both time spent traveling and daily 
path lengths are more dramatic in northern gibbons (Fan & 
Jiang, 2008; Zhang et al., 2014). For example, the daily path 
length of N. concolor in Mt. Wuliang varies monthly from 629 
m to 2 356 m (n=14 months) and daily from 300 m to 3 144 m 
(n=66 days, Fan & Jiang, 2008). The mean daily path length 
of H. tianxing in Mt. Gaoligong varies monthly from 508 m to 1 
926 m (n=14 months) and daily from 345 m to 2 606 m (n=63 
days) (Zhang et al., 2014). Furthermore, travel time can show 
seasonal differences. For example, the H. tianxing group in Mt. 
Gaoligong has been observed to travel (on average) 3.40 h in 
June but 1.27 hin February (Fan et al., 2013). Seasonal impact 
on travel time has also been reported for the two groups of N. 
nasutus in Bangliang (group G1: from 3.92 h in June to 0.52 h 
in March; group G4: 4.92 h in June to 1.37 h in March; Fan et 
al., 2012). 

Both sleeping and resting behaviors in northern gibbons 
are strongly influenced by cold temperatures, with more 
time spent inactive (sleeping and resting) during the 24 h 
cycle in winter (Fan et al., 2012, 2013). In Nankang, Mt. 
Gaoligong, the H. tianxing group spent 20.12 h inactive in 
February, whereas the two groups of N. nasutus in Bangliang 
spent 20.23 h (December) and 19.69 h (January) inactive. 


Gibbons maintain an energy-conserving posture while sleeping 
or resting (sitting with hands clasping their legs: Bicca-Marques 
& Calegaro-Marques, 1998). Northern gibbons usually sit 
huddled together on tree tops in the sunshine during cold winter 
days, and even huddle and share sleeping sites on cold nights 
(Fan et al., unpublished data). 


SOCIAL SYSTEM OF NORTHERN GIBBONS 


Gibbons usually live in small social groups with one adult pair 
and 2-3 offspring (Leighton, 1987). However, some northern 
gibbons live in stable polygynous groups with two breeding 
females (N. hainanus: Zhou et al., 2008; N. concolor: Fan & 
Jiang, 2010a; Guan et al., 2013; Huang et al., 2013; N. nasutus: 
Fan et al., 2010, 2015). Gibbon groups with more than one 
breeding female have also been reported in other species (H. 
pileatus: Srikosamatara & Brockelman, 1987; Hoolock hoolock: 
Ahsan, 1995; H. lar: Sommer & Reichard, 2000), however, 
none of these groups have remained stable for more than two 
years. 

Gibbon groups with more than one adult female were first 
reported in N. hainanus in Delacour’s book On the Indochinese 


gibbon (Xu et al., 1983), and later in N. concolor from Mt. 


Wuliang by Haimoff et al. (1986, 1987). Sheeran (1993) did 
not observe any multi-adult female groups in Mt. Wuliang or 
Mt. Ailao; based on brief observations in the same areas, 
Bleisch & Chen (1991) reported that some N. concolor groups 
included two females who duetted with the adult male, though 
they argued that this could be a sub-adult daughter singing 
with her mother. For a clear demonstration of polygyny, 
close observation of groups with two or more females (with 
young sired by the same male) is required (Bleisch & Chen, 
1991). Jiang et al. (1994a, 1999) suggested that the groups 
they observed were polygynous because the two females 
presented with offspring of the same body size. They thus 
proposed several models to explain the evolution of polygyny 
in N. concolor. First, the large home range occupied by N. 
concolor and N. hainanus might provide sufficient space for 
more individuals in one group. Second, mutual tolerance 
of females or “weak territoriality in females” might facilitate 
polygyny. Third, phylogeny might play a role in the evolution 
of polygyny in Nomascus gibbons. They also argued that 
N. concolor has a relatively large semi-pendulous scrotum 
(Groves, 1972), which should promote higher sperm production 


and, consequently, support a higher level of sexual activity. 


However, the above findings raise questions of why and how 
these gibbons occupy a large home range and females tolerate 
each other. These questions remained unanswered at the 
time as no gibbon groups had been habituated for detailed 
behavioral observation. 

In 2006, Fan et al. (2006) reported that all five groups in a 
small subpopulation of N. concolor at Dazhaizi, Mt. Wuliang, 
included two adult females who gave birth or carried infants 
during the two-year research period. However, females in 
the same group bred at different times, perhaps to reduce 
reproductive competition. During the study, the researchers 
habituated one group and were the first to document the social 


behavior between two females living in the same group. They 
reported that the two females fed in the same tree, sang 
synchronously, groomed each other, and groomed and were 
groomed by the adult male, but never slept in the same place 
(Fan et al., 2006). They even shared meat when one female 
killed a flying squirrel (Fan & Jiang, 2009) and cooperated 
in evicting a “floating” (unmated) female from their territory 
(Fan et al., 2006). Further longitudinal observations showed 
that the two females repeatedly bred with the resident group 
male (Huang et al., 2013) and maintained similar spatial 
relationships to him (Fan & Jiang, 2010a). More recently, Guan 
et al. (2013) demonstrated that these groups consisted of two 
breeding females even after replacement of breeding males 
and females, and females actively cooperated in maintaining 
their social relationship, rather than co-existing merely through 
tolerance or avoidance. 

The rediscovery of a N. nasutus population has provided 
further opportunity to study the evolution of polygyny in gibbons. 
This population lives in a degraded karst forest patch (Fan 
et al., 2011a), which is distinct from the habitat of any other 
gibbon population. Most groups in this population include two 
adult females. Fan et al. (2015) documented that N. nasutus 
benefits from living in larger groups because these groups can 
invade the home range of smaller gibbon groups and access 
more non-fig fruit, though larger groups also spend more time 
traveling. Female cooperation to evict unmated females has 
also been observed in this population but no group consisted 
of more than two breeding females, suggesting that territories 
defended by a male could not support a third female (Fan et al., 
2015). 

Long-term population monitoring can provide insight to 
better understand the sociality of northern gibbons. Although 
floating males in N. concolor populations have been observed 
regularly at Dazhaizi (Fan et al., unpublished data), females 
in the polygynous groups do not disperse and form new 
groups with these floating males. Two newly formed N. 
nasutus groups at Bangliang consisted of two adult females 
(Wei et al., 2017; Shao-Gan Wei, personal communication), 
though solitary males were observed frequently (Fan et al., 
2010, 2015). Similar results have been reported from a 
small population of N. hainanus; although solitary males were 
observed living in the site, one of the two newly formed groups 
consisted of two females (Bryant et al., 2016; Deng et al., 2017). 
These results reject the hypotheses that females are forced to 
live in their natal groups by delayed dispersal or are forced to 
immigrate into monogamous groups because suitable habitat 
and mates are not available (Bleisch & Chen, 1991; Liu et al., 
1989). Consequently, Fan et al. (2015) proposed that polygyny 
is an evolutionarily stable mating system in certain habitats and 
for some species. 

To secure their resource requirements, northern gibbons 
need to maintain a large territory in their heterogeneous 
habitats (Liu et al., 1989; Fan & Jiang, 2008, 2010b; Zhang 
et al., 2014). Within each territory, important food species 
occur in large patches, which can provide adequate food for two 
breeding females and their offspring, but perhaps not enough 
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for three females (Fan & Jiang, 2008; Fan et al., 2015). The 
ability to digest substantial amounts of leaves in response to 
fruit scarcity (Fan et al., 2009a, 2012) enables these gibbons 
to minimize feeding competition within groups when fruit is 
not available and forage in large polygynous groups. Adult 
Nomascus males have a large scrotum and females in the 
same group breed at different times (N. hainanus: Zhou et al., 
2008; N. concolor: Huang et al., 2013), which both reduce 
reproductive competition and enables polygyny (Path | in 
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Figure 2). Southern gibbons live in more homogenous habitats 
(e.g., S. syndactylus: Chivers, 1974; H. lar: Bartlett, 2009; 
Brockelman et al., 2014), with smaller home ranges (Bartlett, 
2011). Although they utilize large fruit trees (Bartlett, 2009), 
intensive interspecific competition depletes food resources, 
which may promote within-group competition (Marshall et al., 
2009) (Path Il in Figure 2). That does not mean food patch size 
in the tropical gibbon habitat is physically smaller than that for 
northern gibbons, but that on average, food availability is less. 
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Males benefit from forming polygynous groups through 
increased mating opportunities (N. hainanus: Zhou et al., 2008; 


N. concolor: Huang et al., 2013; N. nasutus: Fan et al., 2015). 


The cooperative female relationships documented above (N. 
concolor: Fan et al., 2006; Fan & Jiang, 2010a; Guan et al., 
2013; N. nasutus: Fan et al., 2015) suggest that females in 
polygynous groups also benefit from larger group sizes. The 
potential benefits to females include resource access from 
successful intergroup encounters (Fan & Jiang, 2010a; Fan 
et al., 2015). Resident females in the same group appear to 
produce similar-aged offspring. After maturing, these offspring 
can disperse together, which might increase their survival rate 
and probability of establishing their own territory (Fan & Jiang, 
2010a; Fan et al., 2015). In addition, females might benefit 
from improved inclusive fitness if females born within the same 
group disperse together (Fan et al., 2015). Even if females 
do not obtain any direct or indirect benefits from polygyny, 
bi-female groups could form if the cost of evicting the second 
female exceeds the cost of tolerance (Fan & Jiang, 2010a; Fan 
et al., 2015). 

We need to note that H. tianxing living in northern montane 
forest shares many similarities in ecology with the three 
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Nomascus species (reviewed above), but no stable group of 
more than one adult breeding female has been reported. A 
possible reason is that food patch size in the habitat of H. 
tianxing is smaller than that of Nomascus species (Path III in 
Figure 2). Given that stable polygynous groups have only been 
reported from Nomascus species living in different habitats, 
phylogeny may also play an important role in the evolution of 
polygyny in gibbons (Bryant et al., 2015; Di Fiore & Rendall, 
1994; Jiang et al., 1999). 


FUTURE RESEARCH DIRECTIONS 


Standardized methods 

Comparative research using the same or similar methods will 
be extremely useful for understanding gibbon ecology and 
social system evolution. First, the same methods should 
be applied to quantify forest structure, plant abundance and 
distribution, and food availability (Vogel & Janson, 2011). 
Second, the same or similar behavioral definitions and 
methodologies should be used to record behavior. For 
example, several researchers distinguished foraging as “the 
group moved slowly, spread out over a broad front, feeding 
briefly from many of the trees as they passed through them” 


from a feeding bout when all group members fed in the same 
tree at the same time. Kappeler (1981) and Ahsan (2001) 
recorded foraging, but other researchers did not separate 
foraging from feeding (Bartlett, 2009; Fan et al., 2008, 2012, 
2013; McConkey & Brockelman, 2011). Palombit (1997) used 
focal animal observations to record behavior, whereas other 
researchers have used scan sampling (Ahsan, 2001; Chivers, 
1974; Fan et al., 2008, 2012, 2013; Gittins, 1982; Lappan, 
2007a). Even for researchers using scan sampling, scan 
intervals often differ (e.g., 5 min: Bartlett, 2009; Fan et al., 
2008, 2012, 2013; 10 min: Ahsan, 2001; Gittins, 1982; 30 min: 
McConkey & Brockelman, 2011). Furthermore, comparisons 
would be strengthened if researchers used the same method 
to analyze data. Although most researchers present diet and 
time budget results as percentages (Ahsan, 2001; Bartlett, 
2009; Fan et al., 2008, 2009a; Gittins, 1982; McConkey & 
Brockelman, 2011), Fan et al. (2012, 2013) used time devoted 
to different activities and food categories because active time 
showed considerable disparity as a consequence of variations 
in day length at high latitude areas. It would strengthen 
the ability to compare results if all authors published their 
data as both percentages and actual times for both diet and 
time budgets. Finally, altitudinal range must be considered 
when calculating home range size and daily travel length 
because both variables can be impacted at steep study sites 
(Monterroso et al., 2013). 


How do northern gibbons defend their large home range? 
Gibbons are territorial and regularly and aggressively defend 
their home range, which is thought to enable exclusive access 


to a feeding territory (Bartlett, 2009; Brockelman, 2009). 


Northern gibbons occupy home ranges several times larger 
than those of tropical gibbons (N. concolor: Fan & Jiang, 
2008; N. nasutus: Fei et al., 2012). However, they do 
not travel any longer than gibbons that occupy smaller home 
ranges. Gibbons use their loud morning singing to defend their 
territories. However, northern gibbons do not increase singing 
to defend their larger home ranges (Fan et al., 2009b; Yin et 
al., 2016). In fact, hoolock gibbon groups appear to decrease 


their singing in areas with low group density (Yin et al., 2016). 


Therefore, the way in which northern gibbons defend their large 
territories still requires further investigation. 


Social relationship in groups with different social systems 
Intragroup social behavior in gibbons has received little 
attention (Bartlett, 2011). From the limited information 
available, males have been observed grooming females more 
often than the reverse in the monogamous groups of S. 
syndactylus (Chivers, 1974; Palombit, 1996), H. lar (Palombit, 
1996), and H. hoolock (Ahsan, 2001). Guan et al. (2013) 
studied the grooming patterns of two polygynous groups of 
N. concolor and found that females groomed males more 
often than the reverse, a pattern also observed in other 
polygynous primates (Ahsan, 2001). Studying intragroup social 
relationships in populations with different social systems will 
help to better understand the evolution and maintenance of 
social systems in gibbons. 


Are northern gibbons adapted to digesting leaves? 
Northern gibbons eat more leaves than do tropical gibbons (but 
see S. syndactylus in Bartlett, 2011, Table 17.2). However, 
how they resolve the chemical and mechanical challenges 
associated with leaf eating (Lambert, 1998), especially during 
cold winters with reduced or no fruit, remains unclear. Thus, 
comparative research on their digestion system, including their 
microbiome, will help clarify their ecology and social systems. 


Kinship among females 

Female kinship may be a factor affecting social relationships 
(Clutton-Brock & Lukas, 2012). Females living in polygynous 
groups may benefit through inclusive fitness if they are kin, 
and is highly possible due to the previous observation of two 
N. nasutus females dispersing together (Fan et al., 2010). 
However, as it is difficult to observe group formation or female 
displacement in the wild (Brockelman et al., 1998), the use 
of non-invasive genetic methods may help determine kin 
relationships among females (Lappan, 2007b). 


How do northern gibbons find food within their large home 
ranges? 

Do northern gibbons use spatial memory to find food resources 
in the same manner as tropical gibbons (Brockelman et al., 
2014)? In northern habitats, gibbons appear to prefer ripe 
fruit pulp (Fan et al., 2009a, 2013). Scarce but predictable 
food resources led to the development of the “ecological 
intelligence hypothesis” (Milton et al., 1981), that is, primates, 
such as great apes, which rely on ephemeral and scattered 
fruits need larger ranges and have evolved advanced cognitive 
capacities that facilitate optimal food finding strategies. Does 
the “ecological intelligence hypothesis” explain the foraging 
behavior of gibbons? Long-term comparative research across 
gibbon species and between gibbons and large-bodied apes 
will help elucidate northern gibbon foraging behavior in their 
larger home ranges. 


Comparison between foraging behavior 

Given the different interspecies food competition at different 
sites, tree diameter at breast height (DBH) or crown size 
may not be accurate indicators of food availability. We highly 
recommend using feeding party size, feeding bout length, and 
number of feeding patches visited as indices to represent 
food availability, especially as gibbons often deplete their food 
patches. However, this prerequisite needs careful investigation. 
If these indices are validated, then detailed comparison among 
the foraging behaviors of gibbons living in different social 
systems can be made. 


Life history and cultural behavior 

Gibbons born in polygynous groups may develop better social 
skills for maintaining relationships due to their interactions 
with more social partners while growing up. Therefore, 
ontogeny may also play a role in the maintenance of polygyny. 
Female gorillas familiar from immaturity usually exhibit friendly 
relationships in adulthood (Harcourt, 1979). Studying the 
development of social relationships in gibbons is difficult 
because of their long lifespans, and thus longitudinal data is 
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required to answer the above question. Only with long-term 
study and the accumulation of data can we reveal unusual 
observations, such as “dance” (Fan et al., 2016) and “juggling” 
behaviors (Deng & Zhou, 2016). Current evidence implies that 
there is a great deal more to learn about gibbon species, and 
a wide variety of behavior remains to be explored in future 
studies. 
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ABSTRACT 


In primate species with social systems consisting 
of one-male breeding units (OMUs), resident male 
takeover represents a major challenge to individual 
reproductive success and mating strategies. The 
golden snub-nosed monkey (Rhinopithecus roxellana) 
is characterized by large multilevel societies (MLS) 
comprised of several OMUs and all-male units (AMUs); 
however, the factors and mechanisms associated with 
resident male takeover, which offer important insight 
into primate reproduction and social strategies, are 
still poorly understood. Based on 5-year monitoring 
data from a free-ranging herd of golden monkeys 
from the Qinling Mountains in China, we categorized 
three phases of an OMU, that is, a rising phase, 
developing phase, and declining phase. The rising 
and declining phases were unstable periods in which 
male takeover in an OMU might occur. Factors 
causing takeover, such as leader male rank, fighting 
ability, reproduction rate, and affiliation (proximity, 
allogrooming), were analyzed for males and females 
and for different OMUs. Results indicated that the 
new resident male’s fighting ability was lower than 
that of the former resident male in 23 cases. After 
replacement, the rank order of the new resident 
male significantly declined. Females involved in a 
takeover increased their distance from the resident 
male and decreased mating frequency during the three 
months prior to takeover. Females with infants under 
one-year-old had a marked effect on the specific time 
of takeover occurrence. These results suggested that 
female choice was the main factor deciding whether a 
takeover attempt was successful. Furthermore, rather 


266 Science Press 


than male conflict, females more often initiated and 
affected takeover and outcome, implying that the social 
status and competitive ability of the males played 
lesser roles during takeover. 


Keywords: Rhinopithecus roxellana, Female mate 
choice; Takeover; Resident male tenure; Social 
network analysis 


INTRODUCTION 


In mammal species with polygamous breeding groups, a new 
male can enter an established group and forcibly expel the 
resident adult male (Paul et al., 2000; van Schaik, 2000). This 
process is called male replacement or takeover (Rudran, 1973; 
Wheatley, 1982). Resident males monopolize the reproductive 
behavior of the adult females and maintain their long-term 
position in the group after male takeover (Inoue & Takenaka, 
2008). Correspondingly, females may engage in several 
alternative behavioral strategies in response to changes in social 
dynamics during male takeover. For example, female emigration 
can occur immediately after the collapse of a social group when 
the previous male leader is ousted (Pusey & Packer, 1987). 
Alternatively, females may also leave their former unit before 
conflict or takeover, resulting in possible unit collapse. Therefore, 
female transfer and mating choice can play important roles for 
new male mating opportunities (Swedell, 2000). 
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Resident male migration can influence the social organization, 
mating strategies, and genetic structure of a species or 
population (Kuester & Paul, 1999), especially in polygamous 
primates that form multilevel societies (MLS). Only a small 
number of primate species, including hamadryas baboons 
(Papio hamadryas), gelada baboons (Theropithecus gelada), 
and snub-nosed monkeys (Rhinopithecus spp.), are reported 
to live in MLS, within which one-male units (OMUs) travel, feed, 
and rest together to form a cohesive band (Dunbar, 1988; 
Grüeter et al., 2017). In wild hamadryas baboons, most females 
show their acceptance of the intruding male following a takeover 
after a series of behavioral stages (Swedell, 2000). In gelada 
baboons, however, females residing in the same OMU are 
closely related, forming matrilines, so males must aggressively 
compete for reproductive opportunities to maintain their tenure 
as resident leaders (Dunbar, 1988). In the case of snub-nosed 
monkeys, the bands are usually followed by one or more all-male 
units (AMUs) comprised of several adult, sub-adult, and juvenile 
males (Grüeter & van Schaik, 2009). In the Yunnan snub-nosed 
monkey (A. bieti), a resident male’s rank is very important in 
regard to successful takeover (Zhu et al., 2016). 

Golden snub-nosed monkeys (R. roxellana) are an 
endangered primate species, as evaluated by the IUCN, and 
are characterized by a polygamous MLS system (Qi et al., 
2014, 2017). Several OMUs forage and rest together, forming a 
breeding band. Bachelor bands formed by ousted males, both 
former resident males and young males who have reached 
sexual maturity, are usually found nearby. During the mating 
season, the bachelor groups are often observed approaching 
the breeding bands to seek reproductive opportunities (Qi et al., 
2014). As a countermeasure, harem males can form alliances 
to evict solitary males (Qi et al., 2017). Thus, male takeover in 


R. roxellana is influenced by adult males from bachelor groups. 


It has been reported that conflict occurs during takeover in R. 
roxellana (Wang et al., 2004; Zhao et al., 2011). In addition, 
harem female preference can affect the takeover process (Qi et 
al., 2009). However, the specific mechanisms of takeovers, as 
well as the social, reproductive, and demographic factors that 
influence the timing and success of male takeovers, remain 
unclear. Therefore, we carried out a study on the golden 
snub-nosed monkeys from the Qinling Mountains to clarify the 
mechanisms of male takeover and hypothesized that male 
fighting ability had limited impact on takeover events. This study 
was carried by collecting data on male takeover events, social 
ranks, and factors that may affect affiliation between males and 
females before and after takeover, including spatial proximity, 
grooming, reproductive success, and female birth rate. Our 
research provides valuable evidence to help understand the 
mechanism(s) involved in male takeover and reproductive and 
mating strategies in R. roxellana. 


MATERIALS AND METHODS 


Study site 

This research was conducted in the Yuhuangmiao region of 
the Zhouzhi National Nature Reserve (ZNNR), located on the 
northern slopes of the Qinling Mountains (E108°14’—108° 180’, 


N33°45’-33°50’, altitude of 1 400-2 890 m a.s.l.), Shaanxi 
Province, China (Qi et al., 2008). Vegetation in the area 
includes deciduous broadleaf forest (1 400-2 200 m a.s.l.), 
mixed coniferous and deciduous broadleaf forest (2 200-2 600 
m a.s.l.), and coniferous forest (above 2 600 m a.s.l.) (Li & Zhao, 
2007). The region exhibits striking seasonality: the average 
annual temperature is 10.71 °C, with a maximum of 31.51 °C in 
July and minimum of —14.31 °C in January, and average annual 
rainfall is 894 mm (Li et al., 2001). 


Study troop 

Two troops reside in the study area (i.e., East and West 
Range Troops: ERT and WRT) (Li et al., 2000), separated 
by the Nancha River. The WRT consists of the GNG-herd 
and DJF-herd. The GNG-herd was the focal study herd and 
contained a breeding band of 7-13 OMUs and an all-male 
band. Behavioral observations were made at 0.5 to 50 m, 
which allowed behavioral patterns as well as individual ages 
and genders to be recorded under proximity with provisioned 
food. Hand feeding and physical contact between researchers 
and animals were avoided, so all behavioral patterns occurred 
under natural environments and circumstances. 

Individual identification was made according to body 
characteristics, such as pelage coloration, crown hair pattern, 
scars or evidence of previous injury, and shape of granulomatous 
flanges on both sides of the upper lip (Qi et al., 2008, 2009). 


Definition of behavior and takeover 

We defined resident male replacement as behavior involving 
aggressive action, including contacting, chasing, and 
threatening an individual, aggression in response to rejection 
of the attacked, or counter-aggression (Grüeter, 2004). If 
the interval between conflict events between two OMUs was 
longer than 30 s, two conflicts were recorded, otherwise only 
one conflict was recorded. 


Data collection 

The size and composition of the OMUs were recorded in 
March 2008 after all individuals in the herd were identified. 
Other data used in this study were collected from March to 
June, September to the following January, and from August 
2012 to December 2016. In total, data included 36 OMUs 
and 396 monkeys. All individuals in an OMU were scanned 
during the observation period to clarify whether a takeover had 
occurred. Agonistic behavior among OMUs was recorded by all 
occurrence sampling (Altmann, 1974). Allogrooming between 
males and females was recorded. Copulation, involving 
mounting and heterosexual genital contact by intromission and 
pelvic thrusts (Li & Zhao, 2007), was recorded by all occurrence 
sampling, with the initiator and receiver both identified. We 
chose one breeding band as the observation target each 
observation day (1000—1400 h) and allogrooming was recorded 
for 259 d, totaling 906.5 h. 


Data analysis 

Individual social rank was scored according to previous 
research (de Vries, 1998). The birth rate was calculated based 
on the following formula: 
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L 
AEA) , 
ae (1) 
where l; is the number of births observed in the reproductive 
season; F; is the number of adult females during the same 
period; and D is the duration observed and recorded by year 


BR= 


(Eisenberg et al., 1981). 


The reproductive rate in an OMU, as well as male 


reproductive success, was determined as follows: 


OI 
— OF 


(2) 


where / is the sum of infants in the OMU involved in a takeover 
in the first reproduction season; and F is the number of adult 


females without an infant in the replaced OMU. 


The frequency of copulation was determined by C/P/T, where 
C is the time observed during the observation period; T is the 


total observation duration of one OMU, recorded by the hour; 
and nis the number of adult females. 

We used the independent t test to compare differences in 
the reproduction rates of a harem male in different periods of 
tenure. 


RESULTS 


OMU male takeover and male-male conflict 

From August 2012 to December 2016, we recorded 36 
OMUs, of which nine already existed in 2012, 19 underwent 
restructuring during the study period, and eight remained 
unchanged. Of the 36 OMUs, 12 females transferred and four 
OMUs disappeared from the study troop. For the 16 male 
takeover events, only 13 were thoroughly recorded before and 
after, with social rank, tenure, and affiliation between females 
and males documented in detail. Thus, male replacement in 
these OMUs was considered successful (Table 1). 


Table 1 Variation of OMU and transfer of former harem males after replacement 








Form of replacement Replacement Date Rank status Injury Parturition or not* Male transfert 

BZ—TB 2012.8 1618 No a 

FZ LZ 2013.12 10-15 FZ 2/2 0/2 b&a 
PK+>VS 2014.2 213 VS 2/2 0/2 b 

f BB-HT 2013.3 18 No 2/3 1/3 b&c 
Direatreplacement SX>G3 2013.9 1212 3/4 1/4 c 
HB-+BQ 2015.10 a 
XJ—LD 2015.11 No b 
PK-+KO 2016.9 a 

ST-SQ*—HK 2012.12 13814 ST 1/2 1/2 b&c,b&c 
HK-+LZ*—YH 2013.10 14-15 LZ* b,b&c 
Sequential/indirect replacement PK*—BG-+QS 2016.11 

JB>WX-ST = 2013.10-2014.4 9-15 Wx* 1/1 0/1 b 
DXJ-+XC*>G3 2016.11 DXJ* b 


Capital letters in the Replacement column represent the ID of each individuals. *: Resident males got injury during the first replacement. ^: 


Resident male whose tenure was shorter than 10 days. #: The left row was the number of female who was in reproductive season delivery; the 


right row was the female who was in reproduction season did not delivery. +: a, emigrated to other troop; b, transferred into AMU; c, rebuild OMU. 


Resident male takeover accounted for 60% of OMU changes 
in our study. Rank order significantly declined in the new OMU 
compared with the former OMU (independent sample t test: 
t=3.073, df=7, P<0.05), and the reproductive rate of females 
following takeover was not significantly different than that of 
females not involved in a takeover (independent sample f test: 
t=0.244, df=6, P>0.05). 

Of the 13 takeovers, only five involved conflict, totaling 
23 conflict events. Furthermore, the rate of a new resident 
male losing the conflict was significantly higher than of it 
winning (independent sample t test: t=9.349, df=23, P<0.001, 
Figure 1). 


Affiliation of males and females before and after takeover 
With regard to the OMUs in which takeovers occurred, 
affiliation between males and females became weaker in 
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the three months prior to takeover than that observed in 
OMUs assembled for 3-7 years without takeover (independent 
sample t test: allogrooming: t=2.754, df=11, P<0.05; 
proximity: t=5.149, df=11, P<0.001; copulation: f=8.103, 
df=11, P<0.001). Furthermore, for OMUs that experienced a 
takeover and rebuilt within three months, both allogrooming and 
proximity were significantly weaker than that observed in OMUs 
assembled for 3-7 years without takeover (independent sample 
t test: allogrooming: f=3.291, df=11, P<0.001; proximity: 
t=4.052, df=11, P<0.01); however, mating frequency was 
significantly higher (independent sample T test: t=—11.056, 
df=11, P<0.001, Figure 2). 


Resident male tenure affects takeover 
According to differences in the number of adult females and 
infants, as well as the reproductive rates, three distinct OMU 


development periods were defined: (1) rising phase in the first 
year; (2) developing phase from years 1-8, during which adult 
female numbers were stable; and (3) declining phase after 8 
years due to a significant reduction in the number of adult 
females and infants as well as reproductive rates (independent 
sample t test: t=—3.254, df=38; P<0.01). 
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Figure 1 The outcome of replacement fight (five days 
recording after replacement) 

Column A: The replacement fight was not observed; B: the times challenged 
male win during replacement fight; C: the times challenged male lost during 


replacement fight. ***: P<0.001. 
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Figure 2 The affiliation between male and female before and 
after replacement 

Column A: the OMU which undergo the replacement and rebuild within 
three months; B: OMU that build 3-7 years; C: the OMU which before the 


replacement within three months. *: P<0.05. 


From 2001 and 2014, the tenure of four of the resident 
males ranged from 4 to 9 months (average 5.6+2.7, mean+ SD), 
whereas the tenure of the other five was 10 years. Seven OMUs 
stayed in the study troop for more than one year, but less than 
8.5 years (average 3.50+2.43, mean+SD). Takeover frequency 
in the first year (rising phase) and tenth year (declining phase) 
were significantly higher than that in the other years (developing 
period) (Wilcoxon rank sum test: Z=—2.981, P<0.01). 














Months of occurrence 

There was a significant difference in takeover month between 
OMUs containing adult females with infants younger than one 
year and females without young infants (Kruskal-Wallis test: 
x?=11.18, df=3, P<0.05). For OMUs containing females 
with infants under one-year-old, male replacements occurred 
from December to the following March, whereas for OMUs 
comprised of females without young infants, takeovers were 
random throughout the year. 


DISCUSSION 


This study attempted to clarify resident male takeover 
scenarios in golden snub-nosed monkeys. Our study analyzed 
factors involved in resident male takeover based on long-term 
observation and individual identification. 


Limited role of males in takeover 

Our results showed that the new resident male often lost 
the conflict with the former resident male, and social rank 
order of the newly formed unit was significantly lower than 
that of the former unit, suggesting that weak conflict ability 
and lower-ranked males could successfully replace resident 
males in R. roxellana. This result conformed to our previous 
hypothesis. Specifically, males played a lesser role in the 
takeover process as former resident males could win the fight 
but still lose the harem. 

In golden snub-nosed monkeys, resident males in breeding 
bands always face the danger of takeover due to the existence 
of all-male groups. In these groups, males struggle to gain 
breeding opportunities, and thus fighting for takeover occurs 
frequently. However, fierce conflict can lead to substantial 
injury regardless of winning or losing. In the wild, small injuries 
can have severe consequences, even resulting in death, which 
can be a long-term threat to animal populations. Therefore, if 
takeover outcome is decided by fighting ability alone, it would 
be harmful for population development during evolution. Thus, 
R. roxellana males do not display fierce conflict during the 
process of male takeover, as supported by evidence from 
captive populations. Ren et al. (2007) reported that this not 
only reduces the probability of injury, but also lowers the cost 
of male takeover. This behavioral pattern is also seen in other 
species, including wild barbary macaques (Macaca sylvanus) 
and some group-living mammals (Johnstone, 2000; Nonacs & 
Hager, 2011). 


Female choice in takeover initiation 

Females involved in a takeover increased the distance from and 
showed an alienated association with the former resident male 
and decreased mating frequency and mating chances in the 
three months prior to takeover. Thus, our evidence strongly 
indicated that male takeover in R. roxellana was primarily 
associated with female behavior changes. Females have the 
option of direct mate choice by joining or deserting the new 
male, or by leaving when the resident male is the target of 
harassment (Sicotte et al., 2017). Male takeovers are largely 
determined by female transfer and mate choice, which both 
play critical roles in female social and reproductive strategies 
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(Sterck et al., 1997; van Schaik, 1989). Rhinopithecus 
roxellana females expressed their choice by alienating the 
former resident male and mating less often, thereby directly 
initiating and influencing the takeover process. 

We also found that mating frequency was significantly higher 
after the takeover. It may be that females balance the cost 
of wasting investment in having offspring against the risks of 
failing to bond with the new male (Mori & Dunbar, 1985). We 
argue that female mate choice is the primary factor initiating 
male takeover. Given the limited male sexual coercion and 
low levels of female aggression, females appeared to select 
OMUs based on the quality of adult males, which may reduce 
the ability of the resident male to monopolize access and 
increase the opportunity of females to mate with different 
males; however, this requires further research. In our study, 
mating frequency was significantly higher after male takeover 
than that before. This may be a way in which females choose 
mates and increase genetic diversity. Furthermore, mating 
after male takeover helps to establish good sexual relationships 
and assists unit members to accept each other and quickly 
stabilize the OMU. Females likely initiate male takeover to 
reduce inbreeding and ensure that their offspring are sired by 
different males over their reproductive lifetime. 


Resident male tenure affects takeover 

Our results demonstrated that takeover frequencies in the first 
year (rising phase) and tenth year (declining phase) were 
significantly higher than that in the other years (developing 
period). It has been reported that resident male tenure is the 
most important factor that negatively affects male reproductive 
success (Sommer & Rajpurohit, 1989). Females are less likely 
to be inseminated by males with relatively short or long tenures, 
which may lead to potential reproductive competition (Qi et 
al., 2009). Thus, a female more often chooses a male with 
a suitable tenure (Inoue & Takenaka, 2008). In our case, the 
reproductive rate after the eighth year was significantly lower 
than that in the previous years. Relatively short tenure (rising 


phase) showed unstable affiliation and lower reproductive rates. 


Adult females may transfer to a more stable environment by 
takeover to improve female reproduction. Female can also 
lower potential reproductive competition by transferring. 


Strategies of females with infants 

Our results demonstrated significant differences in takeover 
month between OMUs containing adult females with infants 
under one year and females without. For OMUs involving 
females with a young infant, male replacement tended to occur 
from December to the following March, whereas for OMUs 
containing females without young infants, male replacement 
was random throughout the year. 

Females carrying infants less than one year old may also 
influence the seasonal breeding of R. roxellana. Related 
studies show that the golden snub-nosed monkey is a seasonal 
breeding species, with breeding and mating seasons occurring 
in March to May and September to December, respectively 
(Qi et al., 2008). In seasonally breeding populations, the 
optimal time for takeovers is between the breeding and mating 
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seasons (Cords, 2004). However, the timing of resident male 
takeovers in R. roxellana does not support this prediction, as 
also reported for R. bieti (Cui et al., 2006; Kirkpatrick et al., 
1998). However, females with an infant less than one year 
old transferred to new OMUs during the interval between the 
mating and breeding seasons. This may be related to the 
multiple male mating exhibited in this species, which is likely 
a tactic used to confuse paternity as sneaky mating between 
females and non-resident males was observed in the present 
study. Thus, it may be beneficial for females who have a very 
young infant to initiate takeover during this time. 

Although other factors can influence takeover, our results 
indicated that female choice was the main factor triggering 
takeover and outcome. A male’s social status and his 
competitive ability played lesser roles. Female mate choice 
was another driver that decided whether a male had a chance 
at successful reproduction. 
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ABSTRACT 


Play behaviors and signals during playful interactions 
with juvenile conspecifics are important for both the 


social and cognitive development of young animals. 


The social organization of a species can also influence 
juvenile social play. We examined the relationships 
among play behaviors, candidate play signals, and 
play bout termination in Tibetan macaques (Macaca 
thibetana) during juvenile and infant social play to 
characterize the species play style. As Tibetan 
macaques are despotic and live in groups with 
strict linear dominance hierarchies and infrequent 
reconciliation, we predicted that play would be at risk 
of misinterpretation by both the individuals engaged 
in the play bout and by those watching, possibly 
leading to injury of the players. Animals living in 
such societies might need to frequently and clearly 
signal playful intent to play partners and other group 
members to avoid aggressive outcomes. We gathered 
video data on 21 individually-identified juvenile and 
infant macaques (one month to five years of age) 
from the Valley of the Wild Monkeys, Mt. Huangshan, 
China. We used all-occurrence sampling to record 
play behaviors and candidate play signals based 
on an ethogram. We predicted that play groups 
would use multiple candidate play signals in a variety 
of contexts and in association with the number of 
audience members in proximity to the players and play 
bout length. In the 283 playful interactions we scored, 
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juvenile and infant macaques used multiple body and 
facial candidate play signals. Our data showed that 
juvenile and infant Tibetan macaques use a versatile 
repertoire of play behaviors and signals to sustain play. 


Keywords: Social play; Play signaling; Play face; 
Macaca thibetana 


INTRODUCTION 


Play is one of the most conspicuous behaviors in which animals 
engage. Social play combines elements of cooperation, 
communication, and reciprocity in participant actions. Play 
also incorporates behavioral modifications from other social 
contexts, such as agonism, mating, and predation, and 
thus the division between play and non-play is not always 
obvious (Burghardt, 2005). Although play behavior has been 
researched for many decades, it is variously defined depending 
on the field of study from which it is viewed. For example, 
Fagen (1981) defined play as behavior that “functions to 
develop, practice, or maintain physical or cognitive abilities and 
social relationships, including both tactics and strategies, by 
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varying, repeating, and/or recombining already functional 
subsequences of behavior outside their primary context” (p. 
65). Fagen’s definition focuses on play function; however, even 
today the ultimate benefits of play remain controversial. From 
a structural perspective, Burghardt (2005) proposed a working 
definition for identifying play behaviors using five key criteria: 
that is, play behavior (1) has limited immediate function; (2) 
has an endogenous component that is voluntary, rewarding, or 
autotelic; (3) is structurally or temporally different from “serious” 


Table 1 Operational definitions 


Term Definition 


behaviors and is incomplete or exaggerated; (4) is performed 
repeatedly but not stereotypically; and (5) occurs when the 
organism is free of stress or social/physical pressures (i.e., 
the player is in a “relaxed field”). Burghardt (2005) further 
posited that these criteria must be met to label a behavior as 
playful in solitary or social contexts. Fagen’s definition of play 
behavior, coupled with Burghardt’s criteria for identifying these 
behaviors, offer a method to recognize and characterize play 
versus non-play behavior in animals of all ages (Table 1). 





Functions to develop, practice, or maintain physical or cognitive abilities and social relationships, 


Play behavior 


including both tactics and strategies, by varying, repeating, and/or recombining already functional 


sub-sequences of behavior outside their primary context (Fagen, 1981, p65). 





Play signal 


Communicatory behaviors that function to promote, cultivate, and manage social play and 


demonstrate playful intentions (Bekoff, 1974; Fagen, 1981; Yanagi & Berman, 2014a, 2014b). 





Successful bout 


Various solitary and social play behaviors are observed 
across animal species (e.g., turtles, Trionyx  triunguis: 
Burghardt et al., 1996; elephants, Loxodonta Africana: Lee 
& Moss, 2014; domestic dogs, Canis familiaris: Horowitz, 
2009); however, play is particularly important to primates, with 
relevance to social and cognitive development (Martin & Caro, 
1985; Palagi et al., 2007). Therefore, studying play can 
provide important insight on a species’ social relationships, 
how a species develops, and how they relate to cognitive 
abilities. Playing with juvenile and infant (hereafter: juvenile) 
conspecifics is typically the first non-mother activity to occur 
in young animals (Bekoff, 1972; Poirier, 1970). Playing 
with peers gradually increases in frequency, complexity, and 
intensity as juveniles age and their social networks expand. 
The repetitive modification and practice of behaviors (e.g., 
mounting and biting) within the play context may yield both 
short- and long-term benefits in future hunting, mating, or 
social interactions (Bekoff & Allen, 1997). The frequency of 
playful behavior typically declines as juveniles transition into 
adulthood, although adults may maintain playful relationships 
with juveniles (e.g., geladas, Theropithecus gelada: Mancini 
& Palagi, 2009; chimpanzees, Pan troglodytes: Palagi et al., 
2004; Shimada & Sueur, 2014). Less frequently, adult-adult 
play can occur in some primate species in both sexual and 
non-sexual contexts (Pellis & lwaniuk, 1999; T. gelada: Mancini 
& Palagi, 2009; bonobos, P. paniscus: Palagi, 2006; ring-tailed 
lemurs, Lemur catta: Palagi, 2009; Pellis & Iwaniuk, 2000; P 
troglodytes: Yamanashi et al., 2018) but may exhibit variable 
forms and functions in species with different social organization. 
In primates, e.g., Macaca spp., social organization exerts 
a pervasive influence on a variety of behaviors, including 
play (Ciani et al., 2012; Fagen, 1981; Maestripieri, 2004; 


Start of the bout marked by exchange of physical contact, chasing, or other play type or play signal. 


Thierry, 2007). Despotic species, such as Japanese macaques 
(Macaca fuscata), with high rates of aggression, low rates of 
reconciliation, and high levels of nepotism are characterized 
by a competitive, defensive, and low-risk play style (Petit et 
al., 2008; Reinhart et al., 2010). Additionally, mothers in 
despotic species may adopt a more restrictive rearing style, 
and mothers of low-ranking individuals may intervene when 
their offspring attempt to play with the offspring of high-ranking 
females (Maestripieri, 2004). In contrast, species that are more 
egalitarian and characterized by low rates of aggression and 
high rates of affiliation and reconciliation, such as Tonkean 
macaques (Macaca tonkeana), have a more cooperative, close 
contact play style and less restrictive mothers (Maestripieri, 
2004; Petit et al., 2008; Reinhart et al., 2010). These play 
pattern differences appear to vary with social organization and 
exist due to the potential devolvement of play from a friendly 
interaction into an aggressive one that has likely negative 
repercussions for the players. Variability in play frequency 
and form also extends to the individual level, in which an 
individual’s “playfulness” is likely influenced by multiple factors, 
including personality (Lampe et al., 2017; Pellis & McKenna, 
1992), prior experience (Cloutier et al., 2013), opportunity to 
play (Panksepp & Beatty, 1980), and neurochemistry (Siviy et 
al., 2011). However, much of this research is limited to human 
children and laboratory rats; therefore, examining individual 
differences in playfulness across multiple species, including 
primates, is necessary. 

Social play fighting in juvenile animals may influence the 
development of dominance relationships later in life (Pellis 
& Pellis, 1996), and the style of play fighting may also be 
related to a species’ social structure and dominance style 
(Fagen, 1981; Palagi, 2006; Petit et al., 2008). Although gentle 
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play fighting may be used to maintain affiliation, more intense 
rough-and-tumble play may establish a dominance hierarchy 
in post-pubertal juveniles (especially in male-male play bouts) 
by testing the strength of play participants (Palagi et al., 2007; 
Pellis & Pellis, 1996). Thus, play style may be predictive of 
social dominance style, and the balance between cooperation 
and competition within a play bout may be different depending 
on the nature of a species’ social system (Palagi, 2006; Petit 
et al., 2008; Reinhart et al., 2010). An individual may use 
certain behaviors, such as slap or chase, to test, cultivate, or 
stabilize a competitive edge in a play bout (van Leeuwen et 
al., 2011). By testing one’s competitive advantage against a 
peer during play, partners can practice aggressive interactions 
that may be necessary later in life to defend, maintain, or gain 
access to resources. Moreover, the frequency of play fighting 
and use of play signals in certain social systems, such as 
the egalitarian social structure characteristic of adult female 
bonobos (P. paniscus), may indicate the necessity of a flexible 
play style to assess and strengthen social relationships (Palagi, 
2006). 

For many primates, social play is frequently coupled with 
playful signals observed throughout a play bout (Fagen, 1981; 
Yanagi & Berman, 2014a). During play, partners transmit 
and receive signals, which include vocalizations (Biben & 
Symmes, 1986; Kipper & Todt, 2002, Vettin & Todt, 2005), 
body movements, gestures, and/or facial expressions (e.g., 
relaxed open mouth or play face; Bekoff & Allen, 1997; Pellis & 
Pellis, 1996). These signals may function to qualify subsequent 
behaviors as playful, maintain a playful context, and/or help 
to avoid an escalation to aggression, especially when the 
behaviors performed are risky or ambiguous (e.g., play bite, 
play slap, or play fight; Bekoff & Allen, 1997; Burghardt, 
2005; Pellis & Pellis, 1996). Play signals can therefore be 
defined as communicatory behaviors that function to promote, 
cultivate, and manage social play and demonstrate playful 
intent (Bekoff, 1974; Fagen, 1981; Waller & Cherry, 2012; 
Yanagi & Berman, 2014a; Table 1). There is debate that such 
signals are only observed during the context of play and are 
distinct from the behaviors used within a play bout, and also 
predict the occurrence of play (Bekoff, 1974; Fagen, 1981; 
Yanagi & Berman, 2014a). However, many potential signals 
discussed in the literature (e.g., play bow in canids or play 
face in primates) are unclear as indicators of the beginning 
of a play bout; instead, these signals often punctuate the 
bout at different points and are variable in their duration, form, 
and intensity. For example, the play face, also known as the 
relaxed open mouth face, is a frequently-observed play signal 
in primates and is commonly associated with close-quarter 
contact, which occurs during play fights (Palagi et al., 2014; 
Pellis & Pellis, 1996; van Hooff, 1967). The intensity, rate, 
and timing of the play face can change with the intensity or 
behavioral content of the play interaction, thus acting as a 
flexible message to indicate playful intent while the dynamics 
of a bout quickly change (Pellis & Pellis, 1996; Špinka et al., 
2016; Symons, 1978; Waller & Cherry, 2012; Waller & Dunbar, 
2005). This flexibility suggests that the play face is most likely 
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multifunctional and may also work to modulate mood (Pellis et 
al., 2011), invite a third party into the play bout, express playful 
intent to a third party, or act as a reward for both partners for 
playful engagement (Spijkerman et al., 1996). Play signals 
likely perform a critical role in prolonging play duration and 
maintaining multiple players in a bout (polyadic play). However, 
it is possible that certain factors also decrease a signal's 
salience, including increase in bout length, bout intensity, 
misdirection or impairment of the signal, changes in audience 
members, or the addition of more players (Bekoff, 1972). 
Short play bouts may be influenced by the misinterpretation 
of play signals and, similarly, long play bouts with aggressive 
behaviors, such as wrestling, may need to include a higher 
frequency of intention reinforcing play signals (Spijkerman et 
al., 1996). Research reflects the variation and degree of 
frequency, flexibility, and functionality of play signals, which 
also appear to be influenced by species social organization. 
Therefore, possessing a diverse repertoire of signals (such as 
vocalizations and head and body movements), in addition to the 
iconic play face or play bow, may be advantageous for juveniles 
to learn and practice; one might expect multiple play signals to 
have evolved in any given taxon. 

For the genus Macaca (23 species), visual signals during 
play, such as bared teeth or play face, are used differently 
depending on the social organizational grade of the species 
(Scopa & Palagi, 2016; Thierry et al., 2000). For example, 
bared teeth as a play signal might be held longer or 
performed more frequently in one species compared to another. 
Macaca species share similarities in social structure, including 
multi-male and multi-female groups, overlapping home ranges, 
and female philopatry; however, interspecific variation can 
be found in patterns of affiliation, reconciliation, dominance, 
aggression, nepotism, and temperament (Thierry, 1985, 1990; 
Thierry et al., 2000). Macaques are categorized on a 
4-grade tolerance scale, with Grade 1 being highly hierarchical, 
nepotistic, and despotic (e.g., M. mulatta or M. fuscata) and 
Grade 4 being more tolerant or egalitarian (e.g., M. tonkeana 
or M. nigra) (Thierry et al., 2000). Tolerant macaque species 
use play signals interchangeably and redundantly to initiate 
and/or terminate play and self-regulate mood (Pellis et al., 
2011; Scopa & Palagi, 2016). They also engage in rapid 
facial mimicry with their play partner, where one partner 
mimics the facial expression of the other (Scopa & Palagi, 
2016). Such mimicry is an unconscious motor mirror response 
and appears to contribute to increased play duration in more 
tolerant species (Mancini et al., 2013; Scopa & Palagi, 2016; 
see Palagi & Scopa, 2017 for discussion). This flexibility in 
signaling may be related to an unpredictable but cooperative 
and less risky play style. Conversely, in despotic macaque 
species, ambiguous play behaviors and miscommunication 
likely generate increased social repercussions and physical 
risks; therefore, play signals may be particularly inflexible, 
specific, and less interchangeable in such species that exhibit 
high levels of aggression and competitive play style (Scopa & 
Palagi, 2016; Thierry et al., 2000; Yanagi & Berman, 2014b). 
In such cases, play signals may alert both participants and 


third parties that the players are “only playing,” and may be 
performed more often when a third party, such as a related 
adult or juvenile conspecific, is present and participation or 
interference is likely (Pellis et al., 2011). 

In free-ranging despotic juvenile rhesus macaques (M. 
mulatta), play signals predict the imminent occurrence of 
dyadic play (Yanagi & Berman, 2014a). Furthermore, most 
signals are non-randomly associated with various initiations 
of play, indicating that signals are used selectively by the 
players depending upon the play content (Yanagi & Berman, 
2014b). For instance, chase play is often associated with 
a crouch-and-stare signal, and leg-peeks are typically used 
by the play recipient, not the initiator (Yanagi & Berman, 
2014b). For Grade 1 despotic macaques, the specific 
use of signaling might decrease miscommunication and help 
reinforce, clarify, or emphasize the playful intent of the sender, 
whereby indicating that participants are “only playing” alleviates 
potential rising tension within the group. Therefore, one might 
hypothesize that multiple play signals may be used in other 
despotic macaque species to clarify social play while mitigating 


potential associated aggression (Yanagi & Berman, 2014b). 


However, the extent to which these signals are species-specific 
is unknown. 

Tibetan macaques (M. thibetana) are the largest macaque 
species and the most derived of their particular lineage 
(Fooden, 1983; Thierry, 2011). Female macaques reach 
adulthood at approximately four or five years of age, generally 
birth their first offspring between four and six years of age, and 


nurse infants for six to twelve months (Zhao & Deng, 1988). 


In males, adulthood begins at approximately seven years of 
age (Zhao & Deng, 1988). Tibetan macaques reproduce 
seasonally, and the number of offspring sired by a male is 
correlated with his dominance rank (Thierry, 2011; Xia et 
al., 2012). Tibetan macaque social organization consists of 
multi-male, multi-female groups of 15-50 individuals, with a 
female-skewed sex ratio (Berman et al., 2004; Li et al., 2007; 
Thierry et al., 2000; Thierry, 2011). This organization is 
centered on dominance hierarchies and kin-bonded coalitions 
(Thierry, 2011). The dominance rank of females is based on 
matrilines, with a daughter generally attaining the dominance 
rank immediately below her mother but above her older siblings 
(Berman et al., 2004; Thierry, 2011; Zhao, 1997). This 
hierarchy influences intergroup competition among females 
and preferential bonds between kin (Thierry, 2011). Most 
males disperse at adulthood and transfer between groups 
during their lifespan, regardless of dominance rank (Thierry, 
2011; Zhao, 1997). Although group males occupy the top 
ranks, females can occasionally outrank males (Berman et al., 
2004). Adult social relationships influence the socialization 
of immature individuals (Thierry, 2011), with the population 
generally showing a kin bias and linear hierarchies (Berman 
et al., 2004). As such, it is expected that the Grade 2 social 
organization of Tibetan macaques will affect juvenile play via 
third-party adult interference, whereby aggressive or affiliative 
behaviors are used to disrupt and terminate play. However, 
little is known about the play style of Tibetan macaques, 


as there are currently no published studies characterizing 
the specific play behaviors and signals observed during their 
playful interactions. Insight on Tibetan macaque juvenile play 
will provide an important foundation to examine their social 
and cognitive development in relation to social structure and 
organization. 

We characterized play behaviors and identified candidate 
play signals in juvenile Tibetan macaques. Tibetan 
macaques show evidence of despotism, with a Grade 2 
social organization, linear dominance hierarchies, and low 
conciliatory tendencies (Berman et al., 2004). It is expected 
that this despotic social structure will impact juvenile play 
behavior, play bout proximity to adults (i.e., potential third-party 
adult interference), frequency and distribution of play signals 
performed by juveniles, and composition of players within the 
bout (i.e., audience members). Specifically, do juveniles use 
specific candidate play signals in various play contexts based 
on the length of the bout and number of audience members 
in proximity? We therefore tested the following predictions: 
1) due to the despotic dominance style, third-party adult 
interference will end play more often than other forms of play 
bout termination (e.g., withdrawal) (Pellis et al., 2011; Thierry, 
2011; Yanagi & Berman, 2014b); 2) play behaviors and play 
signals will show unequal distribution at the beginning of a 
bout, with play signals occurring more frequently to indicate 
that the subsequent behaviors are playful (Yanagi & Berman, 
2014a); 3) the number of play bout audience members will 
be positively correlated with the number of observed play 
signals due to an increased necessity to convey playful intent 
to more individuals (Bekoff, 1972; Spijkerman et al., 1996; 
Palagi, 2006); and 4) play duration will be positively correlated 
with the number of observed play signals to reinforce, clarify, 
and emphasize playful intent (Pellis & Pellis, 1996; Spijkerman 
et al., 1996; Yanagi & Berman, 2014b). We also analyzed 
the distribution of playful behaviors and signals across all 
individuals and investigated the possible effects of sex and age 
on the frequency of dyadic engagement in playful behaviors 
and play signaling. 


MATERIALS AND METHODS 


Subjects and study site 

Video data were collected on 21  individually-recognized, 
free-ranging infant and juvenile Tibetan macaques of the 
Yulingkeng A1 (YA1) group located at the Valley of the Wild 
Monkeys in the Huangshan Scenic District, Anhui Province, 
China (see Berman & Li, 2002 for more information about the 
site). These individuals were between one month and five 
years old based on Macaca age categorization (Thierry, 2011) 
and current group structure maintained by the researchers of 
Anhui University. At the time of the study (summer 2015), 
10 juvenile males were present, including (4) five year olds, 
(3) three year olds, (1) two year old, and (2) individuals less 
than 12 months old; and 11 juvenile females were present, 
including: (1) five year old; (3) three year olds; (4) two year olds, 
and (3) individuals less than 12 months old. Data collection 
focused on juveniles as adult-adult play is rare and may display 
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different frequencies of play signals than that which occurs 
during immature play (Mancini & Palagi, 2009; Palagi et al., 
2004; Shimada & Sueur, 2014). The members of YA1 have 
been habituated to human presence since 1986 for scientific 
research and since 1992 for tourism (see Berman et al., 2004). 
The macaques at this site are provisioned with corn three to 
four times daily by the park staff, and the feedings are visible to 
tourists (see Berman & Li, 2002). The monkeys occasionally 
interact with people (McCarthy et al., 2009). We obtained 
research approval from the Central Washington University 
Institutional Animal Care and Use Committee (#A021507), 
and our research protocols followed the legal requirements of 
the People’s Republic of China and the American Society of 
Primatologists’ Principles for the Ethical Treatment of Primates. 


Video data collection 

Video data were collected at the study site from 3 August to 
19 September 2015 from approximately 0800 h to 1800 h each 
day, resulting in 48 d in the field and approximately 400 h of 


Table 2 Definitions and components of play behaviors 


Type of play Behavior component 


video data. During preliminary data collection, inter-observer 
reliability in individual identification (x=0.86) and intra-observer 
reliability using play behavior and signal ethograms were 
established (play behavior: «=0.93, play signals: K«=0.88, 
bout termination: «=1.00, actor identification: «x=0.87, and 
audience member identification: «=0.86). All-occurrence 
sampling was used to efficiently record high quantity playful 
interactions (Altmann, 1974). All playful interactions defined by 
Fagen (1981) and outlined by Burghardt (2005) were recorded, 
including play behaviors (see Tables 1 and 2) and candidate 
play signals (see Tables 1 and 3), using a Canon HD Vixia 
camcorder. Observations were undertaken from tourist-viewing 
platforms and near feeding sites and other locations where 
juvenile and adult macaques were visible. If a second play bout 
occurred simultaneously while observing a preceding play bout, 
only play behaviors and signals in the initial bout were recorded. 
After the players disengaged from the play bout, recording 
continued until 10 s had elapsed without bout re-initiation. 


Definition 





Chasing Leaping, running, walking 


Locomotive actions, such as running, climbing, and 
leaping towards or away from another individual, in 
which animals alternate roles of chaser and chasee, 
without having body contact with each other. 





Embracing, holding, hugging, 


Slightly resembles wrestling, but in an extremely 
mild form, i.e., holding each other with very slight 











Cuddling i ; À ; 
touching pushing of body, but without any body displacement. 
Often resembles embracing. 
Play in which animals grapple and place their 
Biting, dragging, embracing, E 3 ae : E , 
i ne , : mouths on each other’s body. Typically involves 
ae grabbing, hitting, leaping, lying, For , , : 
Play biting eae i : i similar behavior patterns to wrestling but occurs with 
pinning, pulling, pushing, rolling, A A oe, é ‘ f 
y A ; . biting. Biting and avoiding being bitten with body 
running, tackling, touching, walking : key 
displacement are central activities. 
EN : , Two animals hit each other with their hands for a 
. Hitting with hands, touching, ý i f 
Slapping i ave period without proceeding to a clearer form of play, 
visual fixation cide 
nor terminating the play encounter. 
Dragging, embracing, grabbing, Also known as rough-and-tumble play. Includes 
Wrestling hitting, leaping, lying, pinning, play behavior patterns in which two monkeys 


pulling, pushing, rolling, running, 
tackling, touching, walking 


engage in mutual grasping, pushing, pulling, 
and rolling, without attempts to bite each other. 





Ethogram for rhesus macaque (Macaca mulatta) play: Yanagi & Berman (2014a; 2014b). 


Video data analysis 

From the video footage, the timestamp, actor identity, audience 
member identity, and play signals and behaviors of the 
participants were coded (QuickTime Player for Mac). A 
modified macaque ethogram (Yanagi & Berman, 2014b) was 
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used to analyze the play bout video data, recording all play 
behavior and candidate play signals observed throughout all 
dyadic and polyadic play bouts. Behaviors previously unlisted 
by Yanagi & Berman (2014b) were added to the play signal 
ethogram and included play threat and slap-and-play face 


(Table 3). We estimated each player’s proximity to group 
members as within arm’s reach (<50 cm) or beyond arm’s 
reach (>50 cm). If a group member (adult or juvenile) was 
in proximity to the play bout then he or she was considered 
part of the play bout’s audience, regardless of their orientation 
or activity. Each member of the play bout was counted to 
generate an audience member tally. If an individual performed 
a play signal or behavior without a group member in proximity, 


Table 3 Definitions of play signals 


Play signal Definition 


the audience member tally was categorized as zero. In this 
study, we considered play bouts “successful” when the start of 
the bout was marked by the exchange of contact, chasing, or 
other play behaviors/signals (recorded as “start of play”) (Table 
1). A successful play bout was considered “terminated” when 
a player engaged in 1) non-play activities, e.g., grooming; 2) 
withdrew from the bout; or 3) adult interference occurred (Table 
4; recorded as “end of play”). 





Crouch-and-stare 


Animal's ventral surface is on/near ground and its limbs are fixed, 


while maintaining visual fixation on partner (Symons, 1978). 





Dangle-and-stare 


Animal stares at partner while hanging from an object by hind limbs (Levy, 1979). 





Gamboling 


Bobbing, high stepping gait in which forequarters and hindquarters are alternately 
raised (Symons, 1978). Often accompanied by rotation of the head (Sade, 1973). 





Animal hides behind an object and then peeks at partner, 
alternating the two behavior patterns. 


Hide-and-peek 





Animal stares at partner through its legs with the top of its head against the ground 








Leg-peek : í 
(Symons, 1978). Animal may hold its ankles or place forearms on ground. 
Look-back Animal's body is oriented away from partner in a fixed position on all fours, while the 
head is turned toward the partner over the shoulder (Levy, 1979; Symons, 1978). 
Play face Relaxed, open mouth face, typically observed during play bouts (Levy, 1979). 





Roll-onto-back 


Animal rolls and lies on its back and stares at partner (Levy, 1979). 


-and-stare 





Play threat 
(candidate signal) 


Animal directs a lunge <2 body lengths towards another individual, 
ending the movement by hitting the ground, without facial expression. 





Slap-and-play face Animal hits another individual’s body while simultaneously 


(candidate signal) 


Adapted from Yanagi & Berman (2014b). 


Table 4 Definitions of play bout termination 


Mode of termination Definition 


directing an open mouth face towards the individual. 





Players begin to engage in any behavior/ 


Non-play activities 


activity not considered under the category 


or criteria of play (Berman et al., 2004). 





Players move out of proximity from each 


Withdraw 


other (out of arms reach), and no subsequent 


play behaviors or signals are observed. 





Play bout is interrupted by an adult group member 


Adult interference 


performing aggressive or non-aggressive behaviors 


towards any player (Berman et al., 2004). 


Adapted from Berman et al. (2004). 
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Statistical analysis 

Using IBM SPSS Statistics (Version 23), Vassar Stats Website 
for Statistical Computations (©Richard Lowry, 1998-2013), and 
UCINET (Borgatti et al., 2002), we tested each prediction 
with «=0.05. We only analyzed successful bouts and the 
signals that were displayed during such bouts. We used 
chi-square goodness of fit tests to assess predictions 1, 2, and 
3. Spearman's rank correlation tests were used to analyze 
the average number and rate (average number/min) of play 
signals across bout length (prediction 4). Additionally, we used 
chi-square goodness of fit tests to examine the distribution of 
the total frequency of play signals across varying audience 
numbers. We also used chi-squared goodness of fit tests to 
analyze the distribution of play behaviors and signals for each 
individual. Spearman’s correlation coefficient was used to test 
the correlation between the number of observed play signals 
across audience member categories. We also used MR-QAP 
matrix-based regression analyses in UCINET to analyze the 
effect of sex and age classes on the dyadic frequency of playful 
behaviors and signals. 

While these analyses might be less robust than others, 
chi-square goodness of fit tests are valuable as they can use 
entire datasets to determine overall patterns (Reinhart et al., 
2010). We used chi-squared analyses to investigate general 
patterns in Tibetan macaque play. However, as other authors 
deviate from investigations of general patterns, more complex 
analyses would be preferable to account for the wide variation 
in assumptions that are violated by play-typical data (e.g., 
GLMM). 


RESULTS 


In total, we recorded and analyzed 283 dyadic (=2 individuals) 
and polyadic (>2 individuals) play bouts, with 136 observations 
of the start of play, 183 observations of the termination of 
play, and 94 observations of a complete play bout (where the 
start and termination of play were clearly marked). In the 94 
completed play bouts observed, average length was 64.7 s, 
with a range of 1 to 585 s. The number of players present 
in a play bout ranged from 1 to 5, and all 22 juveniles were 
observed to engage in play at least once. 


Distribution of play across individuals 

Only positively identified individuals were included in the 
analyses (n=19). During play bouts, 16 juveniles exhibited 
both play behaviors and signals during bouts. We 
used chi-square goodness of fit tests to investigate the 
distribution of individual engagement in play behaviors and 
signals. Play behaviors were differentially (as opposed to 
equally) and significantly distributed across individuals (n=19; 
x2(18)=4 016.06, P=0.001). Each juvenile engaged in playful 
behaviors between 4 and 589 times (125.35+169.58). Play 
signaling was also significantly and differentially (as opposed 
to equally) distributed across individuals (n=19; ¥?(15)=561.52, 
P=0.001). Each juvenile exhibited play signals between 1 and 
101 times (25.56+30.93). These results indicate wide variation 
in the frequency and use of playful behaviors and signals 
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across individuals. 


Distribution of play across age and sex classes 

We used MR-QAP regression analyses in UCINET to 
investigate possible correlations between the dyadic 
frequencies of play behaviors, play signals, age, and sex. 
Being of the same age class had no effect on either the 
dyadic frequency of playful behaviors (Y=—0.007X+0.86, 
P=0.41, r2=0.00) or play signaling (Y=—0.008X+4.53, P=0.43, 
r?<0.001). However, being of the same sex was significantly 
related to both the dyadic frequencies of playful behavior 
(Y=1.22X+0.54, P=0.025, r?=0.023) and play signaling 
(Y=16.91X+2.73, P=0.031, r?=0.017). Although these results 
showed that sex had a significant effect on play behavior and 
signaling, the sizes of these effects were small (r?<0.023). 
Given these results, we did not divide the subsequent analyses 
into different dyad-specific sex or age classes. 


Play bout termination 

We used chi-square goodness of fit tests to compare the 
three play termination categories (Table 4) and test prediction 
1 (third-party adult interference will end play more frequently 
than other forms of termination). Results showed a statistically 
significant difference from the expected values (7?(2)=53.52, 
P=0.001), indicating that termination occurred more frequently 
due to non-playful behaviors (n=74) and withdrawal (n=94) than 
to third-party adult interference (n=16). Therefore, prediction 1 
was not supported. 


Play signals 

In total, 415 play signals in playful interactions were coded. The 
average number of play signals seen in completed play bouts 
(n=94) was 1.6 play signals per bout. We used chi-square 
goodness of fit tests to compare the total frequency of play 
bouts initiated by either a play behavior (Table 2) or a play 
signal (Table 3). Results showed a statistically significant 
difference from the expected values, and thus prediction 2 was 
not supported (y2(1)=45.88, P=0.001): play bouts were most 
frequently initiated by play behaviors (n=108) rather than play 
signals (n=28) (an unequal distribution across bouts). In the 
play bouts observed (n=283), nine candidate play signals were 
recorded, including crouch-and-stare (n=41), dangle-and-stare 
(n=21), gamboling (n=3), hide-and-peek (n=1), look-back (n=1), 
play face (n=263), roll-onto-back-and-stare (n=19), play threat 
(n=17), and slap-and-play face (n=49) (Table 5). Six of the 
seven candidate body signals observed in rhesus macaques 
(Yanagi & Berman, 2014a) were also observed in Tibetan 
macaques (crouch-and-stare, dangle-and-stare, gamboling, 
hide-and-peek, look-back, and roll-onto-back-and-stare). Two 
additional candidate signals were also observed (play threat 
and slap-and-play face). 

We compared the frequency of play signals observed for 
zero (n=14), one (n=137), two (n=159), three (n=86), four 
(n=16), and five audience members (n=3) and predicted that as 
the number of audience members increased, the frequency of 
play signals would also increase (prediction 3; Table 5). Results 
showed a significant statistical difference from the expected 


values (x?(5)=335.46, P=0.001); however, the Spearman's we identified when play signals were used by a current player 


rank correlation between frequency of play signals and number as a new member joined or a member withdrew from a 
of audience members was not significant (r(4)=-0.371, P>0.05; bout. Results showed a significant deviation from the expected 
Figure 1). Thus, these results suggest the trend between values (x2(1)=6.025, P=0.001), and most players did not 
play signals and audience members was not linear, but the use a play signal when a new member joined the bout (no 
frequency of play signals across audience member categories signal=230, signal=26) or when an existing member withdrew 
was distributed non-randomly. To examine these results further, (no signal=101, signal=24). 


Table 5 Play signals observed for each audience member category 


Audience Play Dangle- Crouch- Hide- Look- Roll-onto-back- Play Slap-and- Total 























Gamboling 
number face and-stare and-stare and-peek back and-stare threat play face (n) 
Zero 0 1 8 0 0 0 5 0 0 14 
One 84 11 18 1 1 1 7 7 7 137 
Two 104 6 11 1 0 0 4 7 26 159 
Three 61 3 3 1 0 0 2 1 15 86 
Four 11 0 1 0 0 0 1 2 1 16 
Five 3 0 0 0 0 0 0 0 0 3 
Total (n) 263 21 41 3 1 1 19 17 49 415 
We used Spearman’s rank correlation to test prediction 4, as aggressive interactions. In despotic species such as 
that play duration will be correlated with the rate of observed Tibetan macaques, observers of play may interpret a friendly 
play signals. Results showed a strongly significant positive exchange as an aggressive one if it does not contain clear 
correlation between the observed rate of play signals and bout and frequent play signals, and consequently may terminate or 
duration (Spearman’s rank correlation: r(5)=0.991, P=0.001; disrupt the play bout. However, our results did not support 
Figure 2) and the average number of play signals and bout this prediction. In fact, play termination in juvenile Tibetan 
duration (Spearman’s rank correlation: r(5)=0.964, P=0.001; macaques usually occurred because individuals withdrew from 
Figure 3), thus indicating that as the length of a play bout the bout or engaged in non-play behaviors. We observed 
increased, the average number and rate of play signals also only 16 cases in which other monkeys terminated the playful 
increased. interaction: one case by a low-ranked adult male and one case 
by a low-ranked adult female, two cases by a mid-ranked male 
DISCUSSION and five cases by a mid-ranked adult female, three cases by a 
We examined the behaviors and candidate signals used in high-ranked adult male and two cases by a high-ranked adult 
juvenile Tibetan macaque play bouts to characterize the play female, and two cases by a young adult female. We have 
style of this species. Yanagi & Berman (2014b) found insufficient data to draw conclusions regarding the identities 
that juvenile rhesus macaques use play signals in selective of the individuals who terminated these play bouts, but our 
ways, as signals were non-randomly associated with various anecdotal evidence does not follow our prediction both in terms 
subsequent play behaviors, and further hypothesized that of the overall number of third-party terminations or the identities 
signals may be necessary to reinforce and clarify playful intent. of the individuals who terminated the bouts. It is possible 
Despotic macaques may use play signals more frequently to that juvenile macaques avoid areas occupied by adults to 
emphasize the affiliative, rather than aggressive, nature of their directly manage the end of play themselves and avoid adult 
behaviors. Our study showed that juvenile Tibetan macaques interference (Self et al., 2013). Additionally, although prediction 
used similar play signals as rhesus macaques, but because 1 was not supported, due to their despotic dominance style 
these behaviors did not frequently precede the beginning of and based on the Burghardt (2005) criteria for play, young 
a play bout, their use may reinforce and clarify playful intent macaques may avoid adults to maintain play in a relaxed field. 
rather than signal the initiation of a bout. However, further comparative research is needed to test the 
Berman et al. (2004) hypothesized that due to their despotic hypothesis that social organization in all macaque species 
nature, Tibetan macaque adults would interfere with play bouts (including Grades 2 and 3) may influence play style (but see 
perceived by third parties as aggressive. Many behaviors used Reinhart et al., 2010), and therefore play termination by adult 
in play, such as bite, wrestle, and chase, are also used in interference. 


aggressive contexts and could potentially be misinterpreted 
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We found that more play bouts began with play behaviors 
than play signals (not supporting prediction 2). Therefore, play 
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signaling in juvenile Tibetan macaques may be used similarly 
to signaling in juvenile rhesus macaques, in which play signals 
are used to clarify, reinforce, and prolong a play bout between 
partners (Yanagi & Berman, 2014b). Furthermore, our data 
showed a possible trend between the number of audience 
members and number of play signals generated in a bout, with 
more signals occurring with one or two audience members 
(prediction 3; Table 5; Figure 1). Several conclusions can be 
drawn from these results. It appears that individuals use play 
face (n=263; Supplementary Figure S1) more than other play 
signals (n=152), regardless of audience number. This indicates 
the salient nature of the play face signal when there is a 
receiver present (Pellis & Pellis, 1996; van Hooff, 1967) or when 
a combative play behavior, such as play slap, may escalate to 
aggression (Palagi et al., 2007). This was further supported by 
the lack of play face observed in the zero-audience member 
category (when no other group member was in proximity 
to the initiator). Therefore, play face may be an important 
communicative tool in polyadic play when a player is near the 
sender. This result suggests that play face may be a signal to 
others in proximity to the bout rather than generated as a signal 
to the players themselves (Pellis & Pellis, 1996; van Hooff, 
1967). However, play face may be a general signal showing 
playful intention, as it can be used in a variety of contexts and 
to clarify play when the dynamics of the bout change (Pellis 
& Pellis, 1996). Moreover, the occurrence of signaling in the 
zero-audience member category (n=14) was much lower than 
all other categories (n=401), and the lack of observed play face 
may be an artifact of this sampling. Additionally, the frequent 
occurrence of play face may be an involuntary artifact of the 
sender's enjoyment of the bout rather than a message for 
players (Spijkerman et al., 1996). 

Only three play signals were observed in the zero-audience 
member category: (1) crouch-and-stare, (2) dangle-and-stare, 
and (3) roll-onto-back-and-stare. This may indicate the need to 
use complex body and facial signals to attract players to begin 
a bout, rather than a face-only signal. For example, juvenile 
chimpanzees often use attention-getting gestures when a play 
partner does not see the signal sender’s play face (Tomasello et 
al., 1989). In this way, a play signal that involves the combination 
of two or more signals may be necessary to amplify or reinforce 
the sender’s message. Our results support this speculation, as 
juvenile Tibetan macaques used multiple play signals from their 
diverse behavioral repertoire to indicate their willingness to play 
with conspecifics. However, a larger sample of zero-audience 
member play bouts is needed for further investigation. 

Additionally, across all audience member categories, the 
leg-peek signal (see Table 3) was never observed. Yanagi 
& Berman (2014b) found that in rhesus macaques, the 
leg-peek was associated with play initiated by the receiver 
and may have served as a play invitation. However, this 
signal was not present in Tibetan macaque play, indicating 
that leg-peek may be a species-specific signal. Similarly, two 
candidate play signals were observed in the Tibetan macaque 
repertoire that were not observed in rhesus macaques: play 
threat (Supplementary Figure S2) and slap-and-play face 


(Supplementary Figure S3). Adult Tibetan macaques use 
threat behavior as part of an aggressive interaction, and 
it is characterized by an open mouth gesture and another 
body movement, such as a ground slap, raised eyebrow, or 
lunge (Berman et al., 2004). Therefore, the play threat signal, 
observed during play bouts, may potentially be an incomplete 
adult threat behavior, with the absence of the open mouth 
facial expression. Play threat may enable juvenile macaques 
to practice components of adult behavior within the context 
of play (Martin & Caro, 1985). However, further research is 
needed to confirm that the candidate play signals observed 
in Tibetan macaque juveniles do not occur under any other 
contexts. Although the characterization of Tibetan macaque 
play signaling differs from that of rhesus macaques, there 
does appear to be some cross-species similarity in both 
play behavior and signaling repertoires. The similarities and 
differences in their play style may reflect social organization 
(Palagi, 2006); however, more comparative research is needed. 

Relative to audience members, we observed the 
combination of slap (a play behavior) and play face (a play 
signal) more than other play signals, such as gamboling, 
hide-and-peek, and look-back (Figure 1). It is possible that 
the playful intention of a slap behavior needs to be clarified 
as play because of its association with aggression (Burghardt, 
1999; Pellis & Pellis, 1996). In this way, the slap behavior and 
play face signal would contradict each other and, therefore, 
might increase the likelihood of the receiver understanding the 
situation as non-threatening. 

Although we expected that more play signals would be 
observed in polyadic play bouts (at least one actor and two 
audience members) (Spijkerman et al., 1996), we did not 
find support for this prediction. Our results indicated that two 
audience members (three players in total) included the largest 
number of play signals emitted for this bout composition. The 
play signals observed in play groups larger than three members 
were less frequent, and the number of play signals used within 
larger group compositions declined as more members were 
added. This result implies a possible threshold for the saliency 
of play signals past a certain audience member number. lIt 
is possible that play bouts with more players increase the 
complexity of the bout, making them harder to manage (Bekoff, 
1972). Furthermore, the play signal message may lose 
salience as the complexity and size of the play group increase 
beyond three players. This breakdown of communication 
may be more readily observed in despotic macaques due to 
their strict social organization and therefore more risky and 
uncertain social interactions. It is also reasonable to conclude 
that play bouts with fewer members reinforce affiliative bonds 
between group members, whereas larger play groups may be 
used to assess individuals’ physical strengths (Pellis & Pellis, 
1996). Therefore, signaling may be crucial in a small play group 
to reaffirm the social context: affiliation rather than aggression. 
Our results also showed that play signals were less likely to be 
emitted when new audience members were added to a bout or 
when a player left a bout. This indicates that when play signals 
are used by juvenile Tibetan macaques, they are used generally 


and in various contexts, as they were observed before, during, 
and when withdrawing from play (Pellis et al., 2011). 

Finally, we examined the length of play bouts in relation to 
the number of play signals used (prediction 4; Figure 2; Figure 
3). Spearman’s rank correlation showed a strong statistically 
significant correlation, thus supporting this prediction and 
indicating the likely importance of play signals in sustaining 
a play bout. Previous literature has shown that play bouts tend 
to be longer when players use the play face (Waller & Dunbar, 
2005), playful facial expressions punctuate bouts rather than 
initiate them (Palagi et al., 2007), third parties join the play bout 
more often when the play face signal is used, and long play 
bouts are typically more intense and contain more play faces, 
especially when ambiguous behaviors are involved (such as 
wrestle and gnaw) (Spijkerman et al., 1996). This increased 
signaling may be critical for maintaining playful interaction, as 
additional evidence indicates that short play bouts are often 
marked by a misinterpretation of signals, thus halting play. 
Our study demonstrated that candidate play signals in Tibetan 
macaques may be used in similar ways, with the longest play 
bouts containing the most play signals. However, further study 
is needed to examine the causal relationship between play bout 
length and the frequency and distribution of play signals. 


CONCLUSIONS 


Juvenile Tibetan macaques maintained playful interactions 
using multiple candidate play signals, which combined body 
and facial gestures, to begin play, encourage the continuation of 
play, and end play. Across all audience member numbers, play 
face was the most frequently observed play signal, possibly 
indicating its salient and general nature. Tibetan macaque 
juveniles utilized play behaviors and signals similar to rhesus 
macaques; however, the differences observed in their play style 
and signal variety may be influenced by the dominance style of 
Tibetan macaques. 
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ABSTRACT 

Comparative studies of sympatric species are 
essential for understanding behavioral and ecological 
adaptation as well as the mechanisms that can 


reduce resource competition to allow coexistence. 


François’ langurs (Trachypithecus francoisi) and 
Assamese macaques (Macaca assamensis) are 
sympatric primate species found in the limestone 
seasonal rainforests of Nonggang Nature Reserve, 
southwestern Guangxi, China. To explore their 
different adaptation strategies, we collected data on 
diet using scan sampling at 15-min intervals. Our 
results revealed that Frangois’ langurs showed a more 


flexible diet composition than Assamese macaques. 


Francois’ langurs increased dietary diversity and 
mature leaf consumption in response to seasonal 
scarcity of preferred young leaves and fruits, whereas 
Assamese macaques relied heavily on young bamboo 


leaves (Indocalamus calcicolus) in most months. 


These variations reflect the differences in digestive 
physiology, morphology, and the temporal and spatial 
distribution of food resources. 


Keywords: Francois’ langurs; Assamese macaques; 
Diet composition; Limestone habitat 


INTRODUCTION 


Comparative studies of sympatric species are essential for 
understanding behavioral and ecological adaptations, as 
well as the structure of animal communities (Fleagle et 
al., 1999; Fleagle, 2013). Many studies have revealed 
considerable variation in foraging strategies among sympatric 
primate species, with food type documented as one of the 
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principal determining factors (Hadi et al., 2012; Nadjafzadeh 
& Heymann, 2008; Porter, 2001; Powzyk & Mowry, 2003; 
Singh et al., 2011; Tomblin & Cradford, 1994). For example, 
in the Peleonan forest of Sumatra, Presbytis potenziani feeds 
more selectively on fruits, whereas sympatric Simias concolor 
feeds predominantly on leaves (Hadi et al., 2012). In addition, 
sympatric primates can adopt different foraging strategies in 
response to temporal changes in food resources (Dammhahn 
& Kappeler, 2008; Porter, 2001; Standford, 2006; Thorén et 
al., 2011; Tutin & Fernandez, 1993). For example, gorillas 
(Gorilla gorilla) in Lopé Reserve, Gabon, rely on non-fruit 
foods when succulent fruits are scarce in the dry season, 
whereas sympatric chimpanzees (Pan troglodytes) continue 
to find and consume considerable quantities of fruit (Tutin & 
Fernandez, 1993). Many of these variations reflect differences 
in morphology, digestive physiology, and the temporal and 
spatial distribution of food resources (Hadi et al., 2012; 
Nadjafzadeh & Heymann, 2008; Powzyk & Mowry, 2003), and 
represent ecological niche separation, allowing the coexistence 
of sympatric primates. 

The François’ langur ( Trachypithecus francoisi) and Assamese 
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macaque (Macaca assamensis) are sympatric species that 
reside on the limestone hills of southwestern Guangxi, China 
(Zhang et al., 2002). Despite differences in morphology 
and digestive system (Chivers, 1994), previous studies have 
documented similar food habits between the two species, 
with both found to be predominantly folivorous but with fruit 
consumed when available (Ahsan, 1994; Chalise, 2003; Hu, 
2011; Zhou et al., 2006). Thus, these similarities imply 


the possibility of food competition between the two species. 


To date, however, no comparative studies have reported on 
the foraging strategies of sympatric Frangois’ langurs and 
Assamese macaques, which is important for understanding 
their coexistence mechanisms. 

In this paper, we compared the diet of François’ langurs 
and Assamese macaques living sympatrically in the limestone 
habitats of southwestern Guangxi, China. We aimed to: 
(1) determine how food resources differ between them; (2) 
investigate how they adjust behavior in response to seasonal 
changes in food resources; and (3) explore whether differences 


in foraging strategies can explain their coexistence. 
MATERIALS AND METHODS 


Study sites and subjects 


This study was conducted from October 2005 to September 
2006 in Nonggang Nature Reserve (E106°42’-E107°4’, 
N22°13’-N22°33’, Figure 1), Guangxi Province, China. The 
reserve comprises three areas, Nonggang (5 426 hm?), 
Longhu (1 034 hm?), and Longshan (3 949 hm?), which are 
separated by farmland and villages. The reserve consists of 
limestone hills ranging from 400 m to 600 m a.s.l. (Guangxi 
Forest Department, 1993). The vegetation is characterized by 
limestone seasonal rainforest. Annual precipitation was 1 373 
mm (October 2005—-September 2006) and 952 mm (October 
2006-—September 2007). There are two distinct seasons: a 
rainy season from April to September with >80 mm monthly 
rainfall and a dry season in the remainder of the year with <80 
mm monthly rainfall (Zhou et al., 2007). 
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Figure 1 Map of Nonggang Nature Reserve showing the study site and surrounding area 


Our study site is found within the northwestern portion of 
Nonggang (Figure 1). We selected one group of François’ 
langurs (n=9) and two groups of Assamese macaques (Group 
1, n=15; Group 2, n=12) who ranged nearest to our temporary 
camp. The home ranges of the three study groups overlapped 
with each other. The François’ langur group consisted of one 
adult male, five adult females, and three infants. Assamese 
macaque Group 1 consisted of two adult males, four adult 
females, four adult individuals of unidentified sex, and five 
juveniles, and Group 2 consisted of two adult males, four adult 


females, two adult individuals of unidentified sex, and four 
juveniles. 


Ecological sampling 


We conducted vegetation surveys in the main study area at 
the onset of behavioral data collection. We used a stratified 
random sampling method for the placement of vegetation plots. 
We placed 13 plots (50 mx10 m) in the main study area, 
including four in valley basins and nine on hillsides. The 
plots covered most vegetation types described by Shu et al. 


Zoological Research 39(4): 284—290, 2018 285 


(1988). Within the plots, all trees with >5 cm diameter at breast 
height were tagged. In total, we monitored 312 trees from 
30 families at monthly intervals, and recorded the presence of 
young leaves, fruits, and flowers. The relative abundance of 
different plant parts was expressed as a percentage of trees 
bearing the plant parts of interest each month, regardless of 
the size of the canopy. Huang et al. (2010) described monthly 
phenology changes. 


Data collection 

We conducted behavioral observations of the François’ langur 
group for 126 d (7-22 d each month) and of the Assamese 
macaque groups for 58 d (3—9 d each month). Each day, data 
collection began when the monkeys were first encountered, 
and ended when they disappeared or entered their sleeping 
sites. We used scan sampling with 15-min intervals. Each 
scan lasted 5 min, followed by 10 min of inactivity until the 
next scan began. We recorded the activity of each individual 
seen during each scan. We watched each individual for 5 s 
after detection. The behaviors included four activity categories: 
resting, moving, feeding, and social behavior (e.g., grooming 
and playing). To avoid sampling bias toward certain individuals 
or a particular age-sex class, we collected behavioral records 
on as many different individuals as possible during a scan so 
that all individuals in the focal group were included, but we 
sampled no individual more than once. When the individual 
was feeding, we recorded plant species and parts eaten, e.g., 
young leaf, mature leaf, fruit, flower, seed, or other (e.g., petiole 
and stem). During the study period, a total of 8 168 behavior 
records were obtained from 3 052 scan samples for Frangois’ 
langurs, which included 1 599 feeding records. We collected 6 
525 behavior records from 1 666 scan samples for Assamese 
macaques, which included 1 259 feeding records. 


Data analysis 
Because few records were collected in September or October 
2005 for the Assamese macaque groups, we excluded data in 
these months from later analyses. We also excluded records 
for dependent infants and juveniles because these animals 
were not acting independently. We determined the percentage 
of different plant species in the diet of each study group 
by calculating the percentage of feeding records devoted to 
them among annual total feeding records. We calculated 
the Shannon-Weaver diversity and evenness indices based on 
consumption of plant species. The Shannon-Weaver diversity 
index was calculated using the equation: H=—)~P;x\|nP; 
(where P; is the proportion of feeding records of the ith plant 
species) and the evenness index was calculated using the 
equation: J=H’/Hmax=H'/Inn (where n is the number of species 
eaten). Food category composition was expressed as the 
percentage of different plant parts in the monthly diet of the 
study group using monthly total feeding records. Annual food 
category composition was obtained by averaging the monthly 
percentages. 

We used the Wilcoxon signed-rank test to examine 
inter-specific variations in the overall pattern of use of different 
food categories. The Mann-Whitney U test was used to 
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examine seasonal variations in the use of different food 
categories. We used Spearman rank correlations to test 
the relationship between the abundance and consumption 
of different plant parts. All tests were two-tailed, with a 
significance level of 0.05. 


RESULTS 


During the study period, we identified 92 plant species 
consumed by Francois’ langurs and 69 plant species consumed 
by Assamese macaques. Major foods (accounting for >1% 
of total feeding records) contributed to a large proportion of 
the total diet (Frangois’ langur: 80.6%; Assamese macaque: 
85.2%, Table 1). However, the François’ langurs used more 
plant species as major foods than did the Assamese macaques. 
Annual diversity and evenness of use of major food plants 
were higher for François’ langurs (H=2.563, J=0.817) than for 
Assamese macaques (H=1.164, J=0.468). This indicated that 
Assamese macaques concentrated on fewer food species than 
François’ langurs: the top 10 food species accounted for 82.8% 
of the Assamese macaque diet, but only for 54.3% of the 
François’ langur diet. In particular, Indocalamus calcicolus, a 
small bamboo species, contributed to 62% of the annual diet of 
Assamese macaques. 

There was no significant difference in monthly number 
of plant species eaten by Frangois’ langurs and Assamese 
macaques (Z=—0.356, n=10, P=0.722, Figure 2). However, 
inter-specific differences in seasonal dietary diversity in 
response to food availability were observed. Francois’ langurs 
consumed more plant species in the dry season than in 
the rainy season (U=4.000, n;=6, n2=6, P=0.025). Monthly 
number of plant species eaten by langurs correlated negatively 
with monthly fruit availability (rs=—0.840, n=12, P=0.001). 
For Assamese macaques, there was no significant seasonal 
difference in monthly number of food species (U=9.000, n;=5, 
no=5, P=0.548). 

In general, both species were highly folivorous, with leaves 
accounting for 71% and 77.4% of the overall diet for François’ 
langurs and Assamese macaques, respectively (Figure 3). 
Fruits contributed to only a small proportion of the overall 
diet (17.4% for François’ langurs and 13.4% for Assamese 
macaques). However, there was marked inter-specific variation 
in annual food category composition (Figure 3). Assamese 
macaques ate more young leaves (Z=—2.701, n=10, P=0.007), 
whereas Francois’ langurs consumed more mature leaves 
(Z=—2.666, P=0.008) and other items (e.g., seeds and petioles, 
Z=—2.521, P=0.012). Moreover, inter-specific variation in 
seasonal dietary shift in response to food availability was 
detected. François’ langurs consumed more young leaves 
in the rainy season than in the dry season (U=4.000, n;=6, 
No=6, P=0.025, Figure 4). The consumption of young leaves 
correlated positively with their abundance (rs=0.865, n=12, 
P<0.001). Conversely, the consumption of mature leaves 
(U=0.000, P=0.004) was higher in the dry season than in the 
rainy season. A significant and negative correlation was found 
between the consumption of young leaves and mature leaves 
(fs=—0.685, n=12, P=0.014). In addition, langurs consumed 


seeds only in the dry season (U=6.000, P=0.022), and the 
proportion to monthly diet varied from 6% (March) to 22% 
(January). For Assamese macaques, young leaves contributed 
to the highest proportion of the monthly diet almost year-round 
(Figure 4). In contrast to François’ langurs, Assamese 
macaques consumed more young leaves in the dry season 
than in the rainy season, even though the difference was not 
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statistically significant (U=11.000, n;=5, nə=5, P=0.754). It 
is worth noting that young leaves of /Indocalamus calcicolus 
contributed to the bulk of the macaques’ total diet (Table 1) and 
young leaf consumption in most months (Figure 5). Macaques 
increased the consumption of this item in the dry season 
compared with the rainy season, but the difference was not 
statistically significant (U=8.000, P=0.421). 
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Figure 2 Monthly number of plant species eaten by sympatric François’ langurs and Assamese macaques at Nonggang 
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Figure 3 Annual diet composition of sympatric François’ langurs and Assamese macaques at Nonggang 


DISCUSSION 


In this study, both primates exhibited marked differences 
in their choice of food species. Compared to François’ 
langurs, Assamese macaques concentrated more foraging 
effort on a few food species, especially Indocalamus calcicolus, 
which accounted for 62% of their diet. This small bamboo 
is superabundant and endemic to limestone hills (Liang et al., 
1988). Thus, using more readily available common plant species 


as a food source may be an effective strategy for Assamese 
macaques to survive in limestone habitats. In contrast, François’ 
langurs showed more flexibility in food choice, as reflected in 
their higher dietary diversity and evenness indices. Moreover, 
François’ langurs tended to select less common plant species 
as favored foods. Of the top 10 food species, only two 
(Pithecellobium clypearia, Clausena anisum) belonged to the 
dominant tree species in the habitat (Zhou et al., 2006). 
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Although both primates were highly folivorous, they exhibited 
dietary differences: François’ langurs consumed more mature 
leaves and less young leaves than Assamese macaques. 
These variations may reflect differences in feeding strategies 
in response to changing food availability. In typical colobine 
species diets, young leaves and fruits are often preferred 
foods (Oates, 1994; Yeager & Kool, 2000). When young 
leaves and fruits became scarcer in the dry season, François’ 
langurs significantly increased their consumption of mature 
leaves. Mature leaves are abundant and available, but rich in 
cellulose and poor in nutrition (Richard, 1985). They usually 
serve as fallback foods for primates during periods of preferred 
food scarcity (Marshall & Wrangham, 2007). In contrast, 
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Assamese macaques maintained a high level of young leaves 
in the diet almost year-round, even an increase in young leaf 
consumption in the dry season. This could be related to the 
fact that Assamese macaques relied heavily on young leaves 
of Indocalamus calcicolus, which contributed to the bulk of their 
diet in most months and exhibited a similar tendency to young 
leaf consumption. Young leaves of this bamboo species are 
plentiful throughout the year, and their availability is invariant 
between seasons (personal observation), which may provide a 
long-term staple food resource for the macaques in limestone 
habitats. Thus, the bamboo-dominated diet could be a key 
factor in the differences in fallback foods between Assamese 
macaques and François’ langurs. 
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Table 1 Plant species used as major foods by sympatric Francois’ langurs and Assamese macaques at Nonggang 





Frangois’ langur 


Assamese macaque 








Species Family Part? F%P Species Family Part? F%> 
Ficus nervosa Moraceae YL,ML,F,FR 9.39 Indocalamus calcicolus Bambusoideae YL 62.02 
Pithecellobium clypearia Mimosaceae YL,ML,S 7.82 Ficus nervosa Moraceae YL,ML,F,FR 4.01 

Aristolochia longgangensis Aristolochiaceae YL,ML,ST 5.94 Guihaia argyrata Palmae F,P 3.64 
Canthium dicoccum Rubiaceae YL,ML,F,FR 5.76 Sinosideroxylon pedunculatum Sapotaceae YL,FR 3.55 
Embelia scandens Myrsinaceae YL,ML 5.57 Canthium dicoccum Rubiaceae YL,ML,F,FR 2.28 
Ventilago calyculata Rhamnaceae YL,ML,F,FR 4.57 Burretiodendron hsienmu Tiliaceae YL 1.55 
Ficus glaberrima Moraceae YL,ML,FR 4.26 Croton euryphyllus Euphorbiaceae YL 1.55 
Pueraria thunbergiana Papilionaceae YL,ML 4.13 Ficus microcarpa Moraceae YL,FR 1.46 
Wrightia pubescens Apocynaceae YL,ML,F,S 3.50 Ventilago calyculata Rhamnaceae YL 1.37 
Clausena anisum Rutaceae FR 3.32 Berchemia floribunda Rhamnaceae YLEP 1.37 
Zenia insignis Caesalpiniaceae YL 3.19 Lepionurus sylvestris Sapotaceae YL,ML 1.18 
Pseudostreslus indica Moraceae YL,ML 2.88 Sapium rotundifolium Euphorbiaceae FR 1.18 
Sabia swinhoei Sabiaceae YL,ML 2.57 

Randia spinosa Rubiaceae YL,ML 2.50 

Sinosideroxylon pedunculatum Sapotaceae YL,FR 2.32 

Urobotrya latisquama Opiliaceae YL,ML,F 2.07 

Cuscuta chinensis Convolvulaceae YL,ML 2.00 

Tirpitzia ovoidea Linaceae YL,ML,F 1.75 

Burretiodendron hsienmu Tiliaceae YL 1.69 

Ficus gibbosa Moraceae YL,ML,FR 1.63 

Boniodendron minor Sapindaceae YL,ML 1.50 

Teonongia tonkinensis Moraceae YL,ML 1.25 

Cudrania cochinchinensis Moraceae YL,ML 1.01 


Parts eaten: YL, young leaf; ML, mature leaf; FR, fruit; S, seed; F, flower; P, petiole; ST, stem. °F%: Percentage of total feeding records. 


As colobines, langurs possess digestive systems distinct 
from the cercopithecine macaques. Langurs have an enlarged 
and sacculated forestomach, which serves as the primary 
fermentation chamber (forestomach fermentation) (Chivers, 
1994). In contrast, macaques have an enlarged caecum 
or colon as the primary fermentation chamber (caecocolic 
fermentation) (Chivers, 1994; Hladik, 1978). The former is 
most efficient for mammals that rely on foods high in structural 
carbohydrates, whereas the latter is most efficient for mammals 
that primarily feed on readily digestible foods (Alexander, 1993; 
Lambert, 1998). Thus, Frangois’ langurs can consume mature 
leaves in large quantities as fallback foods in response to 
preferred food scarcity, whereas Assamese macaques depend 
more on young leaves. 

In addition to digestive physiology, dietary variation may 
be related to differences in positional behavior reflected in 
anatomy (Chalise, 2003; Fleagle, 2013; McGraw, 1998). It 
is likely that different postural abilities enable species to 
exploit different food resources. Although sitting is the 
most common feeding posture for both primate species, 
Assamese macaques adopted the stank/forelimb-suspended 
posture more frequently than Francois’ langurs during feeding 
(unpublished data for François’ langurs, Huang et al., 
2015). The stank/forelimb-suspended posture was used most 
frequently when Assamese macaques fed on young leaves of 
Indocalamus calcicolus. Macaques always stand bipedally on 


the ground and grasp the stem with one forelimb while they 
pluck young leaves with the other. 
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Figure 5 Monthly percentage of feeding records devoted 
to young leaves of Indocalamus calcicolus by Assamese 
macaques at Nonggang 
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In summary, sympatric Frangois’ langurs and Assamese 
macaques at Nonggang adopted different foraging strategies in 
response to changes in the temporal and spatial distribution of 
preferred foods, which appears to reduce competition for food 
resources and allows for sympatry. 
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ABSTRACT 


The Tibetan macaque, which is endemic to China, 
is currently listed as a Near Endangered primate 
species by the International Union for Conservation 
of Nature (IUCN)(2017). Short tandem repeats 
(STRs) refer to repetitive elements of genome 
sequence that range in length from 1—6 bp. They 
are found in many organisms and are widely 
applied in population genetic studies. To clarify the 
distribution characteristics of genome-wide STRs and 
understand their variation among Tibetan macaques, 
we conducted a genome-wide survey of STRs with 
next-generation sequencing of five macaque samples. 
A total of 1 077 790 perfect STRs were mined from our 
assembly, with an N50 of 4 966 bp. Mono-nucleotide 
repeats were the most abundant, followed by tetra- 
and di-nucleotide repeats. Analysis of GC content 
and repeats showed consistent results with other 
macaques. Furthermore, using STR analysis software 
(lobSTR), we found that the proportion of base pair 
deletions in the STRs was greater than that of 
insertions in the five Tibetan macaque individuals 
(P<0.05, t-test). We also found a greater number of 
homozygous STRs than heterozygous STRs (P<0.05, 
t-test), with the Emei and Jianyang Tibetan macaques 
showing more heterozygous loci than Huangshan 
Tibetan macaques. The proportion of insertions and 
mean variation of alleles in the Emei and Jianyang 
individuals were slightly higher than those in the 
Huangshan individuals, thus revealing differences 


Science Press 


in STR allele size between the two populations. 
The polymorphic STR loci identified based on 
the reference genome showed good amplification 
efficiency and could be used to study population 
genetics in Tibetan macaques. The neighbor-joining 
tree classified the five macaques into two different 
branches according to their geographical origin, 
indicating high genetic differentiation between the 
Huangshan and Sichuan populations. We elucidated 
the distribution characteristics of STRs in the Tibetan 
macaque genome and provided an effective method 
for screening polymorphic STRs. Our results also 
lay a foundation for future genetic variation studies of 
macaques. 


Keywords: Tibetan macaque (Macaca thibetana) 
genome; Short tandem repeats; Variation analysis; 
Polymorphism; Next-generation sequencing 


INTRODUCTION 


Short tandem repeats (STRs), also known as microsatellites, 
are highly variable repetitive elements with nucleotide motifs 
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of 1-6 bp and are ubiquitous in most eukaryotic organisms 
(Oliveira et al., 2006). Microsatellite mutations are generally 
derived from replication slippage, leading to the insertion or 
deletion of one or several repeat motifs (Levinson & Gutman, 
1987; Schlétterer, 2000; Schlétterer & Tautz, 1992). STRs 
located in both the coding and non-coding areas of DNA can 
affect gene expression and transcription (Chistiakov et al., 
2006; Li et al., 2004), especially STRs in the coding regions. 
For example, certain human neurological disorders are caused 
by trinucleotide repeat expansions in coding regions (La Spada 
et al., 1991; MacDonald et al., 1993). Previous studies 
have also indicated that the mutation processes of STR loci 
can vary, even among closely related species, resulting in 
interspecific differences in STR allele size (Rubinsztein et al., 
1995; Schlétterer, 1998; Webster et al., 2002). STRs are 
frequently applied to genetic mapping, forensic identification, 
and phylogenetic studies as molecular markers due to their 
high abundance, neutrality, high variability, and codominance 
(Schlétterer, 1998). Recent studies have also emphasized the 
role of STRs in the genetic architecture of quantitative human 
traits (Gymrek et al., 2016). In addition, STRs have been used 
as molecular markers for the study of genetic diversity and 
population genetic structure (Qin et al., 2016; Sunnucks, 2000). 
Here, we determined genome-wide STR variation in Tibetan 
macaque populations based on high-throughput sequencing. 
The Tibetan macaque (Macaca thibetana) is one of 23 extant 
species of the genus Macaca (Primates: Cercopithecidae) 
and is endemic to China. Its geographic distribution extends 
from southwest to east regions in China (Jiang et al., 
1996). Jiang et al. (1996) divided M. thibetana into four 
subspecies according to variations in external morphology, 
hair color, cranial features, and geographical distribution: 
namely, Macaca thibetana thibetana, Macaca thibetana pullus, 
Macaca thibetana huangshanensis, and Macaca thibetana 
guizhouensis. The species is an increasingly studied animal 
model in biomedical research due to its ease of domestication, 
large body size, long life span, and physiological characteristics 
analogous to those of humans (Yi et al., 2012). Studies on 
intraocular pressure (Liu et al., 2011) and liver transplantation 
(Ji et al., 2015) are two good application examples. However, 
habitat destruction, deforestation, and illegal poaching have led 
to considerable fragmentation of natural Macaca populations, 
with a resulting decrease in wild Tibetan macaque populations 
observed in recent years (Jiang et al., 2016). Currently, the 
Tibetan macaque is a vulnerable and endangered species 
(Jiang et al., 2016), and is classified as a Category Il 
species under the Chinese Wild Animal Protection Law. Thus, 
studies on its genetic variability and diversity are urgently 
needed to assess genetic differentiation and develop effective 
conservation strategies. To date, previous studies have 
recognized 56 STR markers in Tibetan macaques by genomic 
library construction and cross-species amplification (Jia et al., 
2011, 2012; Li et al., 2014; Yang et al., 2017). Nevertheless, 
the numbers of identified STR markers are insufficient for 
accurate estimation, conservation, and management of Tibetan 
macaques. Most STRs in the Tibetan macaque genome 
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remain unidentified and there is no information available on the 
comparative analysis of whole genome STRs due to the lack of 
genomic data. 

Whole genome profiles of STR variations in non-model 
organisms have not yet been conducted, despite the discovery 
of STRs in the 1980s (Guichoux et al., 2011). This is due 
to the associated expense as well as the inefficient and 
laborious construction of genomic libraries and enrichment of 
clones containing STR motifs (Castagnone-Sereno et al., 2010; 
Powell et al., 1996). Recently, the advent of high-throughput 
sequencing technology has produced considerable genomic 
data, providing the opportunity for genome-wide analysis of 
STR variations. To date, STR marker development and 
polymorphism analyses based on high-throughput sequencing 
data have been performed in humans (Willems et al., 2014), 
cattle (Xu et al., 2016), maize (Qu & Liu, 2013), pig (Liu et al., 
2017a), faba bean (Abuzayed et al., 2017), and Angelica gigas 
Nakai (Gil et al., 2017). However, the distribution and variation 
of genome-wide STRs in the Tibetan macaque remain largely 
unknown. To identify and characterize all potential STR loci 
and clarify the differences in STRs among Tibetan macaque 
individuals, we systematically profiled STR distribution in the 
Tibetan macaque genome and STR variations between Emei, 
Jianyang, and Huangshan individuals based on whole genome 
resequencing data. These profiling results will contribute to our 
understanding of Tibetan macaque species variability. 


MATERIALS AND METHODS 


Genomic samples and high-throughput sequencing 

Five Tibetan macaques, two obtained from Huangshan 
Mountain (Anhui Province, HS) and three captured from Emei 
Mountain (one individual) and Jianyang city (two individuals) 
in Sichuan Province (SC), China, were used for whole 
genome resequencing. Genomic DNA was isolated from 
blood and muscle tissue samples (Supplementary Table S1). 
All procedures were approved by the local Ethics Review 
Board and were in accordance with all relevant national and 
international regulations. We generated 150-bp paired-end 
reads with insertion sizes of approximately 350 bp on the 
Illumina HiSeq X Ten platform. In total, 93.68-104.19 Gb of 
clean data were obtained after filtering out low-quality reads 
using the NGS QC toolkit with a quality score of 20 and 
percentage of read length less than 75% of given quality (Patel 
& Jain, 2012) at a sequencing depth of over 30x in the five 
samples, respectively. The resequencing genome data were 
deposited in the GSA dataset (http://bigd.big.ac.cn/gsa/) under 
accession number CRA000789. 


Genome assembly and characterization of STRs 

The high-quality clean reads of one Emei Mountain macaque 
were assembled using the short reads assembling program 
ABySS 2.0 (Simpson et al., 2009). The resulting short contigs 
(<150 bases) were excluded from the assembly. We then used 
MSDB v2.4 (Du et al., 2013) to search STRs using the same 
settings and parameters as described previously (Liu et al., 
2017b). 


Genotyping of STRs based on resequencing data 
Genotyping of the STRs was performed using lobSTR software 
(version 4.0.6), a rapid and accurate algorithm for STR profiling 
based on whole-genome sequencing data (Gymrek et al., 
2012). The algorithm can avoid gapped alignment and address 
specific noise patterns in STR calling. In this study, a lobSTR 
reference was first constructed based on rhesus macaque 
(GCF_000772875.2) STR data (generated in Liu et al., 2017b) 
using lobstr_index.py script. The resequencing data were then 
aligned to the rhesus macaque reference genome with default 
parameters. The output BAM files were sorted with SAMtools. 
Finally, we analyzed STR allelotypes based on the alignment 
files of the five samples. In total, we obtained 840 356 STRs 
across the five macaque individuals after removing low-quality 
STRs (QUAL<30), which were used for downstream analyses. 
Statistical analyses were performed using SPSS version 19. 


Genetic relatedness analysis of STRs 

Tetranucleotide STR loci from all five samples were employed 
for genetic analysis. Each STR allele different from the 
reference allele size was coded as a binary string of all 
zeros, and the remaining alleles were set to one. Genetic 
distance was calculated using the vegan package in the R 
program (version 3.0.2) (Team, 2000). Finally, we constructed 
a phylogenetic tree based on the neighbor-joining method in 
MEGA 5.2 (Hall, 2013). 


Screening of polymorphic STRs 
Allele numbers of each STR loci were counted based on 
genotyping data, and high-quality, polymorphic tetranucleotide 


Table 1 Assembly results for the M. thibetana genome 


loci were screened using VCFtools by the application of 
“-min-alleles 4” and “--maf 0.1” parameters. Ten loci from 
these STRs were randomly selected to validate efficiency and 
polymorphism in the 16 Tibetan macaque samples by agarose 
gel electrophoresis and capillary electrophoresis. Cervus 
software was employed to calculate the observed number 
of alleles (Na), mean observed heterozygosity (Ho), mean 
expected heterozygosity (He), and polymorphic information 
content (PIC). 


RESULTS 


Sequencing, assembly, and STR overview of the Tibetan 
macaque genome 


We resequenced the genomes of five Tibetan macaques 
from Sichuan and Huangshan populations using Illumina 
sequencing technology, generating over 30-fold sequencing 
depth for each individual. After data filtering, we obtained 
a total of 312 252 249 clean sequencing reads with 150-bp 
and 350-bp insert fragments, corresponding to 93.68-95.45 G 
base pairs with GC content of approximately 45% and over 
94.06% Q20 bases (base quality more than 20) in each sample 
(Supplementary Table S1). In the assembled genome, we 
acquired a total of 1 223 752 contigs with a total length of 
2.66 Gb and N50 of 4 966 bp. The average GC content 
of the genomic sequences was 40.48%. The maximum and 
average lengths of the contigs were 97 085 bp and 2 172.38 bp, 
respectively (Table 1). The length distributions of the contigs 
can be seen in Supplementary Figure S1. 








Item Number 

Total number of sequences examined (n) 1 223 752 
Total size of examined sequences (bp) 2 658 459 556 
Mean length of examined sequences (bp) 2 172.38 

N50 4 966 

N90 1 108 

GC content (%) 40.48% 
Number of sequences containing STRs (n) 521 523 
Number of sequences containing more than one STR (n) 252 041 


Hereafter, we analyzed the distribution of STRs in the 
assembly using MSDB2.4 software. All 1 223 752 contigs 
generated in this study were used to identify potential STRs 
with minimum repeats of 12, 7, 5, 4, 4, and 4 for all motifs from 
mono- to hexa-nucleotide STRs, respectively. As expected, 
perfect STRs were the most abundant category, followed 
by compound STRs, interrupted compound and interrupted 
perfect STRs, and complex STRs (Supplementary Table S2). A 
total of 1 077 790 perfect STRs (accounting for approximately 
0.78% of the assembly) were mined in 521 523 contigs, with 


252 041 sequences containing more than one STR (Table 1). 


The frequency, average length, and GC content of the potential 
STRs were also analyzed. Among the 1 077 790 STRs, 
mononucleotide repeats with the highest frequency of 234.7 


loci/Mb were the most abundant types, accounting for 57.89% 
of the total number of STRs, followed by tetra- (181 344 loci, 
16.83%) and di-nucleotide (165 769 loci, 15.38%) repeats, with 
hexanucleotide (6 347 loci, 0.59%) repeats the least abundant 
(Table 2). Average length of the STR core sequences increased 
with repeat motifs from mono- to hexa-nucleotide STRs, except 
for trinucleotide repeats, with the total average length found 
to be 19.19 bp. Analysis of GC content for the six STR 
types showed that GC content was highest in dinucleotide 
STRs, accounting for 48.75% of dinucleotide STRs in the 
assembly, followed by tri-, hexa- and tetra-nucleotide repeats. 
The lowest GC-content was found in the mononucleotide STRs, 
accounting for 0.48% in the genome (Table 2). 
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Table 2 Number, length, frequency, density, and GC content of perfect STRs 











Total Average Frequency Density GC GC 
Types Total Counts i 
Length (bp) Length (bp) (loci/Mb) (bp/Mb) Length Content (%) 
Mono- 623 930 10 397 456 16.66 234.7 3 911.083 50 200 0.48 
Di- 165 769 3 646 216 22 62.36 1371.552 1 777 603 48.75 
Tri- 64 529 1 216 338 18.85 24.27 457.535 436 005 35.85 
Tetra- 181 344 4 372 964 24.11 68.21 1 644.924 1407970 32.20 
Penta- 35 871 873 535 24.35 13.49 328.587 230 585 26.40 
Hexa- 6 347 172 866 27.24 2.39 65.025 60 271 34.87 
20 679 37 
Total 1 077 790 19.19 405.42 7 778.706 3962 634 19.16 
(0.78%) 
We also found that most repeat motifs were AT-rich, except ai 
for the dinucleotide repeats. With a frequency of 233.81 loci/Mb, ah 
the A/T mononucleotide repeat was the most frequent motif A 
in the STRs, followed by the AC/GT (42.6 loci/Mb), AG/CT Jaso | 
(17.09 loci/Mb), AAAC/GTTT (13.53 loci/Mb), ATTT/AAAT E$ 
(11.76 loci/Mb), AAGG/CCTT (11.7 loci/Mb), AAAG/CTTT S 100 | 
(11.35 loci/Mb), AAC/GTT (8.47 loci/Mb), GTTTT/AAAAC (6.22 E 
loci/Mb), and ATT/AAT (5.28 loci/Mb) motifs. The frequency of al 
the remaining motifs was 36.52 loci/Mb (Figure 1). The number o ; | E „=... eee = | : 


of major repeats ranged from 12 to 40 for mono-, 7 to 27 for di-, 
5 to 20 for tri-, 4 to 15 for tetra-, 4 to 12 for penta-, and 4 to 8 
for hexa-nucleotide repeat, respectively (Figure 2). These loci 
accounted for 97.89% of the total counts of each STR type. 


Alignment and genotyping of STRs 


To explore differences in STR distribution among the Tibetan 
macaques, we used lobSTR software to examine the 
high-throughput sequencing data of the five individuals from the 
Huangshan and Sichuan populations (Supplementary Table 
S1). A total of 1 060 419 STRs were found, with the 
1 294 455 perfect STRs located in the rhesus macaque 
chromosomes used as a reference. After removing low-quality 
STRs (QUAL<30), we achieved a total of 840 356 STRs across 
the five macaques. On average, we obtained genotypes for 
713 685.8 STR loci, with a mean coverage of 5.94-fold per 
individual, of which 354 834.2 STR loci (49.72%) had 5-fold 
greater call coverage (Figure 3A). Further observation showed 
that 526 699 STR loci (62.68%) of the 840 356 STRs called 
in our study were shared by the five macaques and 35 352 
STR loci (4.20%) were unique to one individual. In addition, 63 
124, 87 390, and 127 791 STRs were identified in two, three, 
and four individuals, respectively (Figure 3A). The STR loci 
showed greatest alignment to chr 1, followed by chr 2, chr3, 
chr 7, and chr 5, and finally chr Y (Figure 3B). Comparison 
of allele number of di- and tri-nucleotide STR motifs indicated 
that AT repeats had the highest proportion in dinucleotide 
STR motifs with greater than one allele, followed by AC, AG, 
and CG repeats. For trinucleotide STR motifs with more 
than one allele, AAG, AAT, and AAC were found at higher 
proportions. Conversely, for STR motifs with only one allele (no 
polymorphism), the proportions of CG and CCG were higher 
than those of the other motif types (Figure 3C). 
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Count of STRs (x103) 





Figure 2 Distribution of STR types in the M. thibetana genome 
by repeat time 


Screening and initial validation of polymorphic STRs 

Analysis of the allele number of STRs among the five 
individuals indicated that loci with only one allele were 
underrepresented in all individuals, accounting for 1.41% 
(11 869 loci) of all loci. Among these loci, trinucleotide STRs 
had the highest proportion, accounting for 6.03% (1 726 loci) 
of the total number of trinucleotide STRs; and mononucleotide 
STRs were relatively underrepresented, accounting for 0.54% 


(3 205 loci) of mononucleotide STRs. STR loci with two alleles 
accounted for the highest proportion in di- to hexa-nucleotide 
STRs. The proportion of STRs with two alleles (22.93%, 
31.51%, 55.88%, 52.12%, 58.09%, and 58.15% from mono- 
to hexa-nucleotide STRs, respectively) increased with motif 











length. Loci with more than four alleles (pSTRs) had the 
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Figure 3 Alignment results of STRs 
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highest proportion in mono- (245 673 loci, 41.56%) and 
di-nucleotide loci (52 180 loci, 40.42%), followed by tetra- (13 
307 loci, 18.84%), penta- (2 560 loci, 14.50%), hexa- (456 loci, 
13.85%), and tri-nucleotide STRs (3 771 loci, 13.17%) (Figure 
3D). pSTRs accounted for 37.83% of the total STRs. 
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A: Number of shared and specific STRs among the five M. thibetana samples. B: Distribution of STRs mapped to reference chromosomes. C: Comparison of 


allele number for di- and tri-nucleotide STR motifs. D: Comparison of allele number for the six types of STRs. 


We also characterized 3 325 tetranucleotide STRs that 


showed at least four alleles with a frequency of over 0.1. 


The coordinate information of these loci can be seen in 
Supplementary Table S3. Ten loci from these STRs were 


randomly selected to validate efficiency and polymorphism. 


Among them, the forward primers of six STR loci that were 
amplified successfully in the 16 samples were labeled by 


fluorescence (FAM/HEX). Finally, four to seven alleles were 
detected across each locus of the six loci. The Ho and 
He values ranged from 0.40 to 0.88 (mean=0.55) and 0.58 
to 0.81 (mean=0.69), respectively, and the P/C ranged from 
0.519 to 0.754 (mean=0.625) for each locus (Table 3). Primer 
information of the six markers is shown in Supplementary 
Table S4. 
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Table 3 Number, length, frequency, density, and GC content of 
perfect STRs 


Locus Na WN Ho He PIC 





JROS 4 16 0.875 0.659 0.590 
JROQ 7 16 0.750 0.808 0.754 
JR10 6 15 0.400 0.577 0.524 
JR12 6 15 0.400 0.811 0.751 
JR18 5 15 0.467 0.602 0.519 
JR20 4 15 0.400 0.701 0.613 


Na, observed number of alleles; N, sample number; Ho, mean 
observed heterozygosity; He, mean expected heterozygosity; PIC, 
polymorphic information content.. 


Variation and genetic relatedness analysis based on 
genotyping data of STRs 

To obtain more accurate and reliable results, we focused 
on the 526 699 loci called in all five macaques, including 
354 577 mono-, 96 370 di-, 20 271 tri-, 41 166 tetra-, 12 042 
penta-, and 2 273 hexa-nucleotide loci. We divided these 
loci into four allelotype categories: homozygous reference 
where both alleles are aligned to the reference; heterozygous 
reference where one allele is aligned to the reference; 
homozygous nonreference where both alleles are not mapped 
to the reference but are the same; and heterozygous 
nonreference where both alleles are different and do not 
match the reference. For mononucleotide repeats, as shown 
in Figure 4 and Supplementary Table S5, the proportion 
of homozygous reference/nonreference loci was lower in 
SC than in HS, whereas, the proportion of heterozygous 
reference/nonreference was higher in SC than in HS. The 
distributional pattern of the four allelotype categories in the 
other repeat types was similar to that in the mononucleotide 
repeats (Figure 4 and Supplementary Table S5). These results 
showed that Tibetan macaques in HS had a higher proportion 
of homozygous reference and nonreference loci than that in 
SC for each repeat type. Conversely, SC individuals had more 
heterozygous reference/nonreference loci (Figure 4). Overall, 
there was a far higher number of homozygous STRs than 
heterozygous STRs in all five Tibetan macaques, and the 
difference was significant (P<0.05, t-test). The percentage 
of homozygous loci was positively correlated with the length 
of the repeat motif (Figure 4). Compared with the reference 
alleles, the proportion of base pair deletions in the STRs 
was greater than that of insertions for the Tibetan macaques 
(P<0.05, t-test), and the proportion of insertions and mean 
variations in STR allele size were slightly larger in SC than 
in HS (Supplementary Table S6). On average, these loci 
demonstrated a 3.4 bp difference, with differences greater than 
5 bp accounting for approximately 18.30% in each individual. 
Furthermore, insertions accounted for 30.60% of loci and base 
pair deletions accounted for 49.83% of loci (Figure 5 and 
Supplementary Table S6). 
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To measure genetic relatedness, 41 166 tetranucleotide 
STR loci called in all five macaques were used for 
relatedness analysis. The phylogenetic tree constructed 
by neighbor-joining based on genotyping data of the 
tetranucleotide loci showed that the three individuals from SC 
were clustered onto one branch and the two individuals from 
HS were clustered onto another branch. On the SC branch, 
the two individuals from Jianyang and one individual from EM 
showed little genetic differentiation (Supplementary Figure S2). 
This result was consistent with their geographical distribution. 
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Figure 4 Comparison of the four allelotype categories for 


Percentage of loci (%) 
a 
So 





each repeat type among the five Tibetan macaques 

(Blue) homozygous reference; (red) heterozygous reference; (green) 
homozygous nonreference; (purple) heterozygous nonreference. Columns 
represent the proportion of the four allelotype categories of A011, J021, 
J031, HT-1, and HT-2 from left to right for each repeat type, respectively. 


DISCUSSION 


Genome assembly and STR distribution of M. thibetana 
Recently, with the development of next-generation sequencing 
technology and advancement in bioinformatics analysis, 
considerable research on genetic variation in primates has 
been conducted based on genome-wide sequencing, including 
on humans (1000 Genomes Project Consortium et al., 2015), 
rhesus macaques (Liu et al., 2017b; Xue et al., 2016; Zhong 
et al., 2016), and cynomolgus macaques (Osada et al., 
2015). To date, however, few genome-scale studies have 
been reported in Tibetan macaques, except for the previous 
identification of 11.9 million single nucleotide variants (SNVs) 
(Fan et al., 2014). Sequencing technology now offers an 
unprecedented opportunity to examine STR loci of good 
quality and STR variations among Tibetan macaques. In the 
present study, based on resequencing data, STR variations 
and characterizations were investigated at the genome-level 
in Tibetan macaques. These data provide a novel genomic 
resource for M. thibetana species and enrich the database on 
genetic variants in macaques. 
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Figure 5 Distribution of allele size differences in STRs from the reference for the five M. thibetana individuals 


The identification and development of Tibetan macaque STR 
molecular markers reported in previous studies were based on 
the construction of genomic libraries (Jia et al., 2011, 2012; Li 


et al., 2014), which is both labor intensive and time consuming. 


Here, we employed resequencing data generated on the 
Illumina platform to characterize the distributional pattern of 
STRs. Analysis of the genome showed an average GC 
content of 40.48%, similar to that reported in other macaque 
genomes (Zimin et al., 2014). The N50 and mean length 
of the assembled genome were 4 966 bp and 2 172.38 bp, 
respectively, sufficient for detecting STR loci. We found that the 
proportion of STRs in the assembled genome was slightly lower 
than the 0.83%—0.88% obtained in other macaques (Liu et al., 
2017b). This may be because short contigs of less than 150 
bases were deleted from the assembly to obtain more reliable 


and accurate STR loci, leading to an incomplete genome. 


Mononucleotide STR loci are the most abundant repeat type 
in most organisms (Sharma et al., 2007), as was also found 
in our study. Dinucleotide and tetranucleotide STRs were the 
second most abundant STRs and showed similar frequencies, 
which is in accordance with that observed in other macaques 
(Liu et al., 2017b). Furthermore, analysis of GC content and 
repeat times for each repeat type also showed consistent 
results with previous study (Liu et al., 2017b). These similarities 
are a good explanation for the relatively high success rate 
of cross amplification among macaques (Engelhardt et al., 
2017). We discovered that the predominant STR loci contained 
more A or T bases, except for dinucleotide repeats in which 
the AC motif was the most abundant, which has also been 
observed in humans (Subramanian et al., 2003). This could 
be due to a reported tendency that the poly(A) stretches in the 
genome might be more biased toward AT base pairs during 
STR evolution and GC-rich genome sequences are more 
easily mutated to produce A-rich repeats than AT-rich regions 


(International Human Genome Sequencing Consortium, 2001; 
Subramanian et al., 2003). 


Characterization and polymorphism of STRs based on 
genotyping data 


Within the 840 356 potential STRs identified, 62.68% of loci 
were called in all five individuals and 4.20% of loci were 
only present in one individual. This may be due to the 
low sequencing coverage, besides true differences among 
different individuals. Notably, however, STRs mapped to the 
chromosomes exhibited similar distributional regularity with the 
distribution of STRs observed in other macaque chromosomes 
(Liu et al., 2017b). For instance, most STR loci were aligned 
to chr 1. Several studies have indicated a positive correlation 
between chromosome size and STR number (Liu et al., 2017b; 
Qi et al., 2015). Thus, the distribution of chromosome size 
in M. thibetana may be similar to that of the rhesus macaque. 
In addition, we noted that GC-rich STR motifs showed lower 
polymorphism, which may be due to the correlation between 
GC-rich sequences and functionality (Benjamini & Speed, 2012). 

Analysis of allele number of STR loci showed that only 
1.41% of the loci showed no variation among the five 
individuals. A similar occurrence has been reported in cattle 
genome research (Xu et al., 2016). Here, loci with only one 
allele were found at highest proportion among trinucleotide 
STRs compared with other STRs. Furthermore, pSTRs 
showed the lowest proportion among tri- and hexa-nucleotide 
STRs. These two results showed that variation in tri- and 
hexa-nucleotide STRs was smaller than that in the other four 
STR types. This could be attributed to the fact that these two 
STR types are more frequent in exons (Qi et al., 2016), which 
are related to gene function (La Spada et al., 1991; MacDonald 
et al., 1993). However, pSTRs were overrepresented in 
general, indicating that high polymorphic STRs were relatively 
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abundant in the Tibetan macaque genome. These loci were 
found at the highest proportion in mono- and di-nucleotide 
STRs, followed by tetra-nucleotide STRs; however, mono- and 
di-nucleotide STRs are unstable and more easily produce errors 
in experiments (Morin et al., 2001; Taberlet et al., 1999). We also 
found the number of tetranucleotide STRs was higher than that of 
tri-, penta-, and hexa-nucleotide STRs. Thus, it is more efficient 
to screen polymorphic loci from tetranucleotide STRs. 

From the 3 325 polymorphic tetranucleotide STRs (>4 
alleles) with allele frequencies over 0.1 detected in this study, 
six loci were used for genetic diversity analysis of the 16 
M. thibetana samples. The allele number, Ho, He, and 
PIC results revealed high genetic diversity within Tibetan 
macaques. Recent study reported 12 polymorphic loci in M. 
thibetana screened from 30 STR sites in humans (Yang et 
al., 2017), though the amplification efficiency was lower than 
these STR loci identified in our study. This result indicated 
that the polymorphic STR loci identified in the present study 
can be used for amplification in the M. thibetana genome 
and polymorphisms of STRs were relatively high. The 
combination of STR loci searched in the assembled genome 
and polymorphic STRs identified in this study for isolation of 
STRs was both efficient and time-conserving. 


STR variation among five Tibetan macaque individuals 


Based on the 526 699 STR loci called in all five samples, 
we found Tibetan macaques in SC had more heterozygous 
reference/nonreference loci than those in HS. This indicated 
that genetic diversity in Tibetan macaques was higher in SC 
than in HS. Previous studies based on mitochondrial DNA 
suggest that overall genetic diversity in M. thibetana is high, 
and that genetic variation is higher in the SC population than 


in the HS population (Li et al., 2008; Zhong et al., 2013). 


Our results are consistent with these conclusions. Analysis of 
SNV distributions in the Tibetan macaque revealed that there 
was a higher number of homozygous SNVs than heterozygous 
SNVs (Fan et al., 2014). A similar occurrence was also 
found in the distribution of STRs, namely the majority of 
STRs were homozygous loci in the Tibetan macaque (P<0.05, 
t-test). Compared with the rhesus macaque reference alleles, 
the proportion of base pair deletions in the Tibetan macaque 
was greater than that of insertions (P<0.05, t-test), consistent 
with the small indels genotyped with the Genome Analysis 
Toolkit (GATK) (Fan et al., 2014). The average length of 
loci was 19.19 bp in the present assembly, which is slightly 


shorter than the 20.14 bp for M. mulatta (Liu et al., 2017b). 


These results indicate that a large proportion of STR loci may 
have longer repeats in rhesus macaques than orthologous 
loci in M. thibetana. This has also been reported from 
human-chimpanzee genomic sequence alignments (Webster 
et al., 2002). The STR loci in each Tibetan macaque showed 


an average difference of 3.40 bp from the rhesus macaque. 


The proportion of insertions and mean variations in STR 
allele size were slightly larger in the Emei and Jianyang 
individuals than in the Huangshan individuals, which may 
reveal differences in allele size of STRs between the two 
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populations. Further studies are needed to examine the effect 
of these potential STR variations. 


Genetic relationship analyses of microsatellites 

The neighbor-joining tree constructed based on tetranucleotide 
STRs classified the five individuals into two different 
branches according to their geographic origin. Little genetic 
differentiation was observed between the Jianyang and Emei 
individuals. Tibetan macaques from Sichuan and Huangshan 
were classified into two different subspecies (M. thibetana 
thibetana and M. thibetana huangshanensis) based on their 
external morphological and anatomical variations (Jiang et al., 
1996), and the neighbor-joining tree confirmed the genetic 
differentiation between the two populations. Previous study 
based on mitochondrial DNA has demonstrated that the two 
populations exhibit significant genetic differentiation (Zhong 
et al., 2013), which may be the result of geographical 
barriers. Specifically, gene flow between the two populations 
is obstructed, leading to genetic differentiation (Zhong et al., 
2013). Owing to the limited number of samples, our results 
provide an initial understanding of the genetic variation of 
STRs in Tibetan macaques, and future studies are needed to 
investigate population genetic variations using more samples 
from distinct populations. 


CONCLUSIONS 


In summary, we performed STR characterization in the Tibetan 
macaque genome and provided a genome-wide atlas of 
microsatellite distribution with next-generation sequencing data. 
The polymorphic STR loci identified from the reference genome 
exhibited good amplification efficiency and can be well used 
for population genetics. This profiling will be conducive for 
future genome-wide genetic analyses of the Tibetan macaque 
at the population scale. Analysis of genome-wide STRs also 
preliminarily revealed variation in the genetic diversity among 
Tibetan macaque populations. This study contributes to our 
further recognition of the genetic variation among Tibetan 
macaque species. 
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